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Abstract We show that optimal L 2 -convergence in the finite element method on 
quasi-uniform meshes can be achieved if, for some so > 1/2, the boundary value 
problem has the mapping property H~ 1+s —» H 1+s for s £ [0,so]. The lack of 
full elliptic regularity in the dual problem has to be compensated by additional 
regularity of the exact solution. Furthermore, we analyze for a Dirichlet problem 
the approximation of the normal derivative on the boundary without convexity 
assumption on the domain. We show that (up to logarithmic factors) the optimal 
rate is obtained. 

Keywords L 2 a priori bounds • duality argument • reentrant corners 


1 Introduction 

The finite element method (FEM) is a widely used numerical technique for ap¬ 
proximating solutions of boundary value problems. It is based on approximating 
the solution by piecewise polynomials of degree k. In the classical case of second 
order elliptic equations with an H 1 -coercive bilinear form, the method is of optimal 
convergence order in the /I 1 -norm. An important tool for the convergence analysis 
in other norms such as the L 2 -norm are duality arguments (“Nitsche trick”). The 
textbook procedure for optimal order convergence in L 2 is to exploit full ellip¬ 
tic regularity for the dual problem. Conversely, this procedure suggests a loss of 
the optimal convergence rate in L 2 if hf 2 -regularity fails to hold. This occurs, for 
example, in polygonal domains with reentrant corners. 

Nevertheless, it is possible to recover the optimal convergence rate in L 2 , if the 
exact solution has additional regularity to compensate for the lack of full regularity 
of the dual problem. More precisely: In this note, we consider a setting where an 
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elliptic shift theorem holds for both the dual and the bidual problem in the range 
[-1,-1 + so] for some so G (1/2,1] (see Assumption 11.11) and show that if the 
solution is in the Sobolev space if fe + 1 +( 1_s o) ; then the extra regularity 1 — so can 
be exploited to recover the optimal convergence rate in L 2 (up to a logarithmic 
factor in the lowest order case k = 1). 

In the second part of this note, we consider the convergence in L 2 of the normal 
derivative on the boundary. We show that the optimal rate 0(h k ) (up to a logarith¬ 
mic factor in the lowest order case) can be achieved, if the solution is sufficiently 
smooth. The proof is based on a local error analysis of the FEM as discussed, 
e.g., in I22U23I- Here, we extract error bounds for the flux on the boundary from 
an optimal FEM estimate on a strip of width 0(h ) near the boundary. Although 
we present the convergence of the flux for an H 1 -conforming discretization, the 
techniques are applicable to mixed methods, [IT], FEM-BEM coupling, [16], and 
mortar and DG methods, [181125] . In fact, the results of the present work lead to 
a sharpening of [ 18] , where convexity of the domain was assumed to avoid the 
analysis of a suitable additional dual problem. The techniques employed here are 
in part similar to those developed in [18] . Nevertheless, they are also significantly 
different since we have opted to forego the direct use of anisotropic norms and 
instead rely on weighted Sobolev norms and the embedding result of Lemma 12.11 

The analysis of the optimal convergence of fluxes has attracted some attention 
recently. Besides our own contributions mum, we mention the works [2ll3lll3| 
where similar results have been obtained by different methods. 


1.1 Notation 

For bounded Lipschitz domains 12 C IR d with boundary r := <912, we employ 
standard notation for Sobolev spaces iT s (12), [Tll21|. We will formulate certain 
regularity results in terms of Besov space: for s > 0, s N, and g e [1, oo] the 
Besov space -Bf g (12) is defined by interpolation (the “real” method, also known as 
A-method as described, e.g., in |211l22| l as 

Bi q (n) = (a L sJ (12), a m (rt)) 0i „ e = s - Lsj. 

Integer order Besov spaces SJ, 9 (12) with n G N are also defined by interpolation: 

Bl q {Q) = (A ra_1 (12), H rl+1 (f2)) 1 /2 iq , n G N. 

To give some indication of the relevance of the second parameter q in the definition 
of the Besov spaces, we recall the following (continuous) embeddings: 

A s+e (12 ) C B|,i(i?) C A s (12) C B|, 0 o(12) C H s ~ e (f2) Ve > 0. 

Of importance will be the distance function Sp and the regularized distance func¬ 
tion Sp given by 

5p(x) := dist(a:, F), Sp(x) := h + dist(a:, F). (1.1) 

Later on, the parameter h > 0 will be the mesh size of the quasi-uniform trian¬ 
gulation. Also of importance will be neighborhoods Sp> of the boundary <912 given 
by 

S D ■■= {x G f?\ Sr(x) < D}, 
with particular emphasis on the case D = 0(h ). 


( 1 . 2 ) 
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1.2 Model problem 

We let 17 C R d , d G {2,3}, be a bounded Lipschitz domain with a polygo¬ 
nal/polyhedral boundary and let on be our model problem: 

— V • (A(i)Vu) = / in 17, u = 0 on 917. (1.3) 

We assume that A and / are sufficiently smooth. Moreover A is pointwise sym¬ 
metric positive definite, and A(x) > «o I for some ao > 0 and all x G 17. As usual, 
E3 is understood in a weak sense, i.e., for a right-hand side / G (/7(}(l7)) , the 
boundary value problem (11.311 reads: Find u G F/q( 17) such that 

a(u,v) := f AVu • Vu = (/, v) Vw G-Ho (17). (1.4) 

J a 

We denote by T : {Hq{Q))' —> i7o(l7) the solution operator. We emphasize that 
the choice of boundary conditions (here: homogeneous Dirichlet boundary con¬ 
ditions) is not essential for our purposes. Essential, however, is the following as¬ 
sumption: 

Assumption 1.1 There exists so G (1/2,1] such that the solution operator f H > 
Tf for | TMp satisfies 

ll r /l|jf 1 +* 0 (O) < C 'll/ll(/f 0 1 ““ 0 (f 2 ))' - C \\f\\L 2 (n)- 

Here and in the following 0 < c, C < oo denote generic constants that do not 
depend on the mesh-size but possibly depend on so. We also use < to abbreviate 
< C. 

Remark 1.2 The present problem is symmetric. As a consequence certain dual 
problems that will be needed below coincide with the primal problem. This will 
simplify the presentation but is not essential. Inspection of the procedure below 
shows that we need Assumption ll.il for the dual problem and the bidual problem 
with weighted right-hand side. ■ 

Let T be an affine simplicial quasi-uniform triangulation of 17 with mesh size h 
and Vh '■= So A {T) C Hq( 17) the continuous space of piecewise polynomials of 
degree k. This space has the following well-known properties: 

(i) Existence of a quasi-local approximation operator: The Scott-Zhang operator 
/£ : H 1 ^) ->• S k,1 (T) of |25] satisfies: 

-If u G H^{Q) then l£u G V h = S kA {T). 

— I k is quasi-local and stable: ||Vz. 2 (jc) ^ where ujk is 

the patch of elements sharing a node with K. 

— if has approximation properties: 

\\V j (u-I^u)\\ L , {K) <h l+1 ~ j \\V l+1 u\\ LHuK) , j G {0,1}, 0 <l<k. 

(1.5) 

(ii) For every v G fl Hq{Q) there holds 

z fl h ll« - z \\m(n) < h 1/2 l|v|| B 3/2 (n) . 

(This follows from property (i) and an interpolation argument using the K- 
method). 
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(iii) The space Vh satisfies standard elementwise inverse estimates: for integer 
0 < j < m < k 

\v\ H ™(K) < Ch~ (m ~ j) \v\ H j( K) MveV h . (1.6) 

The Galerkin method for 1H is then: Find Uh G Vh such that 

a(u h ,v) = (f,v) Vu G 14. (1.7) 


Remark 1.3 The restriction to simplicial triangulations is not essential. Our pri¬ 
mary motivation for this restriction is that in this case the space Vh is known to 
have the above approximation properties, the inverse estimates, and moreover it 
has the “superapproximation property” that underlies the local error analysis as 
presented in [241 Sec. 5.4). ■ 


2 Regularity 

2.1 Preliminaries 

A key mechanism in our arguments that will allow us to exploit additional regu¬ 
larity of a function is the following embedding theorem. 

Lemma 2.1 The following estimates hold, if Q C is a bounded Lipschitz do¬ 
main and z sufficiently regular. 

\\^r 1 ^" +£ z \\l 2 (S7) < \\S r 1 ^ 2+S z\\ L 2^q^ < C e \\z\\ H i/a-€(n), £ G (0,1/2], (2.1) 


\\5- r 1/2 z\\ L , {n) < C|ln/i| 1 / 2 ||^|| B v 2(fi) , (2.2) 

\\Sf 1/2 - e z\\ LH a) < C e h- e \\z\\ Bl , hny e > 0, (2.3) 

\\z\\ LHSh )<Ch 1/2 \\z\\ B ^ (ny h> 0, (2.4) 

\\z\\L 2 {r) < C\\z\\ B i/^ n y (2.5) 


Proof The estimate involving Sr in (12.111 can be found, e.g., in [101 Thm. 1.4.4. 3] 
and (12.411 is shown in [131 Lemma 2.1]. The estimates (12.211 . (12.31) . (12.511 follow from 
ID Sobolev embedding theorems for L°° and locally flattening the boundary T 
in the same way as it is done in the proof of m Lemma 2.1]. For example, for 
(12.511 we note that a local flattening of the boundary T and the ID embedding 

IMll°°(o,i) ^ IM|i a (o,i)IM|jr 1 (o 1 i) imply IMlW) im¬ 

plies the estimate ||z||i,a(.r') < ||z|| B i/ 2 ( f2 ) by |2!l Lemma 25.3]. We recall that for 

half-spaces, the upper bound (12.511 can be directly found in [221 Thm. 2.9.3], see 
also the comment in [33 Sec. 32, Eq. (32.8)]. □ 

One of several applications of Lemma f2. II is that it allows us to transform negative 
norms into weighted L 2 -estimates: 
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Lemma 2.2 For e G (0,1/2] and sufficiently regular z there holds 

\\44 < C e \\6$ +1/2 -‘z\\ L * ia) , -1 < p < 0, (2.6) 

11^11(5^(0))' < C| ln/i| 1/2 ||?r 1/2 ^|U 2( o). (2.7) 

Proof Firstly, we show (1231) : 


II# 


r 2 ll(H 1 /2-e(0))' 


{SpZ, v) 

SUP -77—77- 

veH 1 / 2 ~ e (r?) ll v llif 1 / 2 -®(f 2 ) 


sup 

<jEH 1 / 2 ~ e (f2) 


(5p~ 1/2 ~ e z,Sf 1/2+6 v) 
lkl|jf 1 /2-s(fi) 


<C £ \\5 0 r +1/2 - 6 z\\ LHn) , 


where, in the last step, we employed m of Lemma f2.ll Secondly, (EZD follows by 
the same type of arguments, where the application of no is replaced with that 
of dUD. □ 


2.2 Regularity 

We recall the following variant of interior regularity of elliptic problems: 

Lemma 2.3 Let Q he a bounded Lipschitz domain and z G H 1+l3 (17), /3 G (0,1], 
solve 

-V • (AVz) = / into. 

Then, for a constant C > 0 depending only on ||A|| C0 ,i(T 7 )> a o, P, and Q 

ll# -/, v 2 z|U a(n) < c (K-PfWvw + \\z\\ Hl+ „ w ) . 

Proof The upper bound follows from local interior regularity for elliptic problems 
(see |19t Lemma 5.7.2] or [9] Thm. 8.8]) and a Besicovitch covering argument, see, 
e.g., [61 Section 1.5.2] and |15l Chapter 5]. We refer also to |12[ Lemma A.3] where 
a closely related result is worked out in detail. □ 

2.2.1 Refined regularity for polygons and polyhedra 

It is worth pointing out that neither the structure of the boundary r nor the 
kind of boundary conditions play a role in Lemma 12.31 One possible interpreta¬ 
tion of Lemma 12.31 is that z could lose the Tf 2 -regularity anywhere near T. For 
certain boundary conditions such as homogeneous Dirichlet conditions and piece- 
wise smooth geometries r the solution fails to be in H 2 only near the points of 
nonsmoothness of the geometry. With methods similar to those of Lemma [2.31 one 
can show the following, stronger result: 

Lemma 2.4 Let J2 be a (curvilinear) polygon in 2D or a (curvilinear) polyhedron 
in 3D. Denote by £ the set of all vertices of D in 2D and the set of all edges of 
in 3D. Let 5g be the distance from £. Let z G H 1+l3 (17), f3 G (0,1], solve 11.3\) . 
Then, for a constant C depending only on olq, || A|| c0 ,i( 72 )' an>/ ^> 

II# ,3 V 2 z|| i ,2( 1 7) < Cf3 (j|# ^/||l 2 (C?) + |M|.ff-i+0(fi)) • 
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Proof Follows from local considerations as in Lemma [2.31 The novel aspect is the 
behavior near the boundary away from the vertices (in 2D) and the edges (in 
3D). This is achieved with an adapted covering theorem of the type described in 
Theorems 1X31 1X31 The key feature of these coverings is that they allow us to 
reduce the considerations to balls B = B r (x ) and stretched balls B = B^i +e y{x) 
(with fixed e > 0) with r ~ dist(x,<5) and the following properties: either x G 12 
with B r (x ) C 12 or x G .T and B PI 12 is a half-ball. Local elliptic regularity 
assertions can then be employed for each ball B. □ 

Lemma 1X3 assumes that a loss of H 2 -regularity occurs at any point of non¬ 
smoothness of r. However, the set of “singular” vertices or edges can be further 
reduced. For example, in 2D for A = Id, it is well-known that only the vertices 
of 12 with interior angle greater than 7r lead to a loss of full iL 2 -regularity. It will 
therefore be useful to introduce the closed set M s of boundary points associated 
with a loss of LI 2 -regularity. Before introducing this set, we point out that this set 
is a subset of the vertices and edges: 

Definition 2.5 (iL 2 -regular part and singular part of the boundary) Let 

12 be a polygon (in 2D) or a polyhedron (in 3D) with vertices A and edges £. 

1. A vertex A G A of 12 is said to be H 2 -regular, if there is a ball B e (A) of radius 
e > 0 such that the solution u of lil.3\) satisfies u\s c (A)nn G ii (12) whenever 
f G L 2 (12) together with the a priori estimate ||u|| H 2 (B c (A)nO) ^ C||/||z, 2 (r2)- 

2. In 3D, an edge e G £ of 12 with endpoints A\, A 2 is said to be H 2 -regular if the 

following condition is satisfied: There is c > 0 such that for the neighborhood 
S = ( x ) °f the edge e we have the regularity assertion 

M lsnt 2 £ H 2 for the solution u of m whenever f G L 2 (12) together with the 
a priori estimate |M|ff 2 (sn«) < C\\f\\ L 2 ^ n) . 

Denote by A r C A the set of H 2 -regular vertices and by £ r C £ the set of H 2 - 
regular edges. Correspondingly, let i= A \ A r and £ s := £ \ £ r be the set of 
vertices and edges, respectively, associated with a loss of H"-regularity. Define the 
singular set M s as 

M s := A s {j£ s C P (2.8) 

With the notion of the singular set in hand, we can formulate the following regu¬ 
larity result: 

Lemma 2.6 Let 12 be a polygon or a polyhedron. Let M s be the singular set as 
defined in Definition \2.5\ Then the following is true for any solution z G Hq(I 2) 
of 2/ z G 22 1+/3 (12) for some /3 G (0,1], then with 5 m b '■= dist(-,M s ), there 

holds for some C > 0 depending only on ao, ||A|| c0 ,i( 7 jp /3, and 12, 

II'Sm/v^IU 2 ^) ^ c 0 (|I <5 m//IU 2 («) + • 

Proof The proof is based on local considerations as in Lemma l2.4l We recall that 
not all vertices and edges (in 3D) are included in the singular set M s . This is ac¬ 
counted for by a further refinement of the covering employed. We restrict ourselves 
to the 3D situation. Using finite coverings provided by Theorem IA.6I one may re¬ 
strict the attention to balls B r = B r (x) and stretched balls B = B^ 1+e y{x) (with 
fixed e > 0) with r ~ dist(x, M s ) where one of the following additional properties 
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is satisfied: a) x G I? with B C I?; b) x G Ar and B fi I? is a solid angle; c) x G U£ r 
and B (~l 17 is a dihedral angle; d) x lies in the interior of a face and B n 17 is a 
half-ball. We emphasize that we do not need to consider balls B r (x) with x G A s 
or x G £ s since the covering provided by Theorem IA.6I is such that for every such 
x there is a neighborhood U x of x that is covered by (countably many) balls whose 
radii tend to 0 as their centers approach x. 

□ 


2.2.2 Shift theorems for locally supported right-hand sides 

We have the following continuity results for the solution operator T for our model 
problem (11.31) : 

Lemma 2.7 Let Assumvtion \l.l\ be valid. Then T : {Ho{fI))' —> Bq(17) satisfies 


W T fhizw ^ (Btfwr ( 2 - 9 ) 

\\TfWHWM < CeWS^fh^n), 0 < e < so — 1/2. (2.10) 

In particular, if f G L 2 (17) with supp / C Sh, then 

I|77IIb 2 3/»(„) < Ch 1/2 \\f\\ LHO) , (2.11) 

\\Tf\\ H wi(2) < Ceh 1/2 - s \\f\\ LHn) , 0 < s < so — 1/2. (2.12) 


Proof We follow the arguments of [18l Lemma 5.2). The starting point for the 
proof of (12.91) is that interpolation and Assumption 11.11 yield with 0 G (0,1) 

T : (i7^(l7)) , )fl,oc -»• (H 1+s °(L>), H 1 (n)) e<00 = Bl+^~ e) {f2). 

Next, we recognize as in [181 Lemma 5.2] (cf. [221 Thm. 1.11.2] or [211 Lemma 41.3]) 

((B o 1 - s »(f2)) , ,(B o 1 (f2)) , ) 0lOO = ((H 1 0 ~ s ’>(n),HUn)) eA y 

d ((tf 1 - 80 ^),.ff 1 ^))^)' = (Bl- So( 1 - e \n)y. 

Setting 6 = 1 — l/(2so) G (0,1/2], we get (-E> 2 i i S ° ( ' 1 ^))' = {^ 2^,1 (^)Y and 

B 2 + cY>° < ' 1 ~ e \^) = B?!oo (^7)- The assertion (12.101) follows from the Assumption ll.il 
and (12.61) with /3 = 0. For the bound (12.111) . we argue as in the proof of Lemma f2.2l 
and use (l2Al) . see also [18l Lemma 5.2]. Finally, the proof of (12.121) follows from 
(12.101) and the assumed support properties of /. □ 

We will also require mapping properties of the solution operator T in weighted 
spaces: 

Lemma 2.8 Let Assumption \1.1\ be valid. Then for v G L 2 (f7) 

r(«F 1 «)ll B 3/. (n) < c^nhnsy^vh^a), (2.13) 

||T(5 r 1 n)|| ii -3/2+e(j7) < Cle/i E \\5 r 1/2 v\\ L 2 {n) , e G (0,so — 1/2], (2.14) 
\\T(Sy 1+2s v)\\ H3/2+Hn) < C e \\Sy 1/2+e v\\ L 2 {n) , £ G (0,a o - 1/2]- (2.15) 
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Proof The results follow by combining Lemmas 12.21 and 12.71 □ 

For the analysis of the FEM error on the neighborhood Sh, we need a refined 
version of interior regularity for elliptic problems. The following result is very 
similar to [181 Lemma 5.4] and closely related to Lemma 12.31 

Lemma 2.9 Let z solve the equation 

—V • (AVz) = v in ft. 

Then there exist Ca (depending only on ao, II A|| c o,i( 7 j) and 12) and C\ > 0 (de¬ 
pending only on 12) such that for z G Bi)^(f2), we have 

II^V^IM^) < C A [\/fhi7ij'||z||+ || v / ^r«l|L 2 (r2\s clh )] • (2-16) 

If the right-hand side v satisfies additionally suppu C Sh and furthermore z = Tv, 
then there are constants Ca (depending only on oto, II Allc 0 ' 1 ^)’ and 12) and c > 1, 
c > c' > 1 (depending only on 12) such that for all sufficiently small h > 0: 

(i) If z e B^(I2) then \\&)( 2 \2 2 z\\ L ^ {n \s- oh ) < C A V\ lnh \\\z\\ B 3 / ^ {ny 

(ii) For every a > 0 there holds 

||<5rV 3 2:|| L 2( r2 \ S . h )<CA[||5r lv “ ;2 lli 2 (f2\S' c ^) + ll A llc 1 . 1 (fi)||^rVx||i,2( fi \ So , ji )j. 

(hi) If z G H 3 2 2+S (f2) for some e G (0,1/2), then for some Ca,c > 0 (depending 
on ao, || 12, and e) there holds 

l|V 2 z|| L 2 (fASsh) < C A ,eh- 1,2+£ \\z\\ m „ +c{n) . 

(iv) If Assumption [77?] is valid, then for some C > 0 depending only on the 
solution operator T and Ca,c °f 00) we have || < C||v||. 

Proof of \2.16\) . (0), $v%): [18] Lemma 5.4] is formulated for —A. However, the 
essential property of the differential operator A that is required is just interior 
regularity. Hence, the result also stands for the present, more general elliptic op¬ 
erator (with the appropriate modifications due to the fact that the coefficient A 
is allowed to be non-constant). In the interest of generality, we have also tracked 
in (12.1611 the dependence on the right-hand side v, which was not done in p~8l 
Lemma 5.4]. A full proof can be found in Appendix 0 

Proof of (ml) : This follows again by local considerations similar to those em¬ 
ployed in the proof of [181 Lemma 5.4] and the obvious bound 5r > h on 12 \ Sch- 
A full proof can be found in Appendix 0 

Proof of Hv 1) : In view of CUD, we have to estimate ||~||_H3/2+ 5 ( 1 7p By the support 
properties of v, the bound (12.121) yields ||z||H 3 / 2 +'(fi) < Ch 1 C~ e \\v\\L^(n)- Inserting 
this in (Imll gives the result. □ 
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3 FEM L 2 -e rror analysis 

Let Uh be the FEM approximation and denote by e = u — Uh the FEM error. The 
standard workhorse is the Galerkin orthogonality 


a(e,v) = a(u — uh,v) = 0 Vt 6 Vt. 
We start with a weighted L 2 -error: 


(3.1) 


Lemma 3.1 Let Assumption \T~]] be valid. Assume that a function z G Hq(L 2) 
satisfies the Galerkin orthogonality 


a(z, v) = 0 Vv G Vh- 


Then 


\\S r 1/2+S z\\ LHn) <C s h 1/2+£ \\z\\ HHn) , £ G (0, s 0 1/2] j (3.2) 

¥r 1/2 z\\L*{n) < Ch 1/2 \\nh\ 1/2 \\z\\ H i(ny (3.3) 

Proof The proof follows standard lines. Define if = T{5~f 1+2e z), which solves 


(v,6 


-l+2e 


z) = a(v,ip) Vv G Hq(L2). 


Then we have by Galerkin orthogonality for arbitrary lip G Vh. 

\\5p 1/2+e z\\ 2 L 2 (n) = a(z, ip) = a(z, ip — lip) < C\\z\\ m (a)\\ip - Iip\\ H ^n)- 

From (12.151 in Lemma 12.81 we have \\ip\\H 3 / 2 + E (n) ^ C e \\Sp 1 ^ 2+£ z\\Li(n) so that 
with the approximation properties of Vh we get 


inf H-I^\\h^o) < C £ h N o r 

iipeVh 


1/2+e || Sf 1/2+s z 


I U 2 (t2)- 


This shows ©2). For ([33D, we proceed similarly using the regularity assertion 
(12.131 and the approximation property of Vh■ □ 

Corollary 3.2 Let Assumvtion ll. l\ be valid and the solution u be in H s (f2). s > 1. 
Then the FEM error e = u — Uh satisfies for e G (0, So — 1/2] 


11-5. 


“ 1/2+£ e|| L 2 (n) < C £ h »~ 1/2+e | 


u ||P := min{s, k + 1}. 


The following Theorem 13.31 shows that the optimal rate of the L 2 -convergence of 
the FEM can be achieved also for non-convex geometries if the solution has some 
additional regularity: 

Theorem 3.3 Let Assumption 1 1. I\ be valid. Let the exact solution u satisfy the 
extra regularity u G H k+2 ~ S °{T2). Then the FEM error u — uh satisfies 


l|w - Uh\\ L *(n) < h k+ 1 \\u\\ Hk +2-s 0{n) . 

More generally, if u G s G [l,k + 2 — so], then 

\\u~ Uh\\ L ^o) < /i s-1+So |Mlff»(«), 1 < s < k + 2 - s 0 . 


(3.4) 


(3.5) 
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Proof of Ij3.4\ ): We proceed along a standard duality argument. To begin with, 
we note that the case so = 1 is classical so that we may assume so < 1 for the 
remainder of the proof. Set e := s o — 1/2 G (0,1/2) by our assumption 1/2 < so < 
1. Let w = Te and let Wh £ V/ be its Galerkin approximation. Quasi-optimality 
and the use of (EH) give us the following energy error estimate: 

Ik - Wh\\w (12) £ Jnf Ik - V\ lifi(fi) < /i 1/2+£ |kllH 3 / 2 +'(i2) 

< k /2+£ || e || (Jf i/ 2 -. (n) ), < h 1/2+ °\\e\\ L ^o). (3.6) 

The Galerkin orthogonalities satisfied by e and w — Wh and a weighted Cauchy- 
Schwarz inequality yield for the Scott-Zhang interpolant I k u 

\\e\\h ( n) = a(e,w) = a(e,w - Wh) = a(u - I k u,w - Wh) (3.7) 

< C\\5f 1/2+£ V(u- Iu)\\ LHn) \\Sr /2 ~ £ V(w - w h )\\ L 2 {ny (3.8) 

We get by a covering argument and (12.61) of Lemma 12.11 

\\S~ 1/2+£ \/(u-ljfu)\\ L 2 {n) < h k \\5~ 1/2+e V k+1 u\\ L 2 (s7) 

<h k \\V k+1 u\\ m/2 - H n)- (3-9) 

It should also be noted at this point that in (13. 91) . the weight 2+e can be 
considered as constant in each element K. For the contribution ||<5]/ 2-e V(w — 
Wh )|| L 2 (n) i n (13.8[> . we have to analyze the Galerkin error w — Wh in more detail, 
which will be done with the techniques from the local error analysis of the FEM. 
We split 17 into S c h U(f?\S c h) where c > 0 will be selected sufficiently large below. 
For fixed c > 0, the L 2 -norm on S c h can easily be bounded with (13.61) by 

||^ /2 - £ V(«; - w h )\\ LHSah) < h 1/2 ~ e ||V( W - w h )\\ L 2 {n) < h\\e\\ L 2 (n) . (3.10) 

The term ||5 r ' "V(w — Wh)\\L 2 (fz\s ch ) requires more care. Obviously, S r '~ < 
4 X jl 2 ~ e on 17 \ S c h- We have to employ the tools from the local error analysis in 
FEM. The Galerkin orthogonality satisfied by w — Wh allows us to use techniques 
as described in [241 Sec. 5.3], which yields the following estimate for arbitrary balls 
B r C B r : with the same center (implicitly, r' > r + 0(h) is assumed in (13.111) 1 

II V(u> - W h )\\ L 2 {Br) < ||V(w - IhW)\\ L 2(B rJ ) + r , l _ r \\ w - w h\\L*(B r ,)- (3.11) 

By a covering argument (which requires r' — r ~ c5p(x), where x is the center 
of the ball B r , and c is sufficiently small) these local estimates can be combined 
into a global estimate of the following form, where for sufficiently small ci > 0 (ci 
depends only on 17 and the shape regularity of the triangulation but is independent 
of h ): 

||4 /2 ~ £ V(wt^ Wh)\\ L *(n\s ch ) ^ (3-12) 

j| S 1 r /2 ~ £ V(w - I k w)\\ L 2 (s7 \ Sccih) + ||^“ 1/2_e k - Wfc)|U 2 (rAs colh ). 

This estimate implicitly assumed c\ch > 2 h, i.e., at least two layers of elements 
separate r from 17 \ S ClC h- We now fix c > 2/ci. The first term in (13.121) can 
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easily be bounded by standard approximation properties of Lemma 12.31 and 
Assumption 11.11 

II4 /2 “ e v( w - /£i0)|U a(flNSoaih ) < fc||4 /a - e v 2 H| ia(0) 

<h p| S\! 2 e e.\\ L 2 ^Q) + ||w||h-3/ 2 +'(i 7)] < ^l|e||z, 2 (fl)- 

In the last step, we have to deal with the term || 5~f 1 ^ 2 ~ £ (w — W/i)||l 2 ( 12 \s cc h ) of 
(13.121) . Lemma l3.1l and (13.61) imply 

II 5~ 1/2 ~ £ (w - w h )\\ L *(n\s ecih ) ^ h~ 2z || 5~ 1/2+£ (w - w h )\\ L 2 {n) 

< h~ 2e h 1/2+£ \\w - w h \\ HHn ) < ft||e|| i2( „). (3.13) 

Here we have used the quasi-optimality of the Galerkin approximation with respect 
to the H 1 -norm. 

Proof of 13.5]) : The above arguments show that the regularity of u enters in the 
bound (tT9l) . For u G H 1 { 12), the stability properties of the Scott-Zhang operator 
ifc show 

||^ 1/2+e V(u - I k h u)\\ma) < /i- 1/2+e ||u||«i ( fl). (3.14) 

Hence, a standard interpolation argument that combines (13.91) and (13.141) yields 
\\5p 1/2+E V{u- IhU)\\ L 2 {n) < h~ 1/2+s+s ~ 1 \\u\\ H s (a) for s e [1, fe + 2 - s 0 ]- Com¬ 
bining this estimate with the above control of ||<5^/ 2-e V (w — u>h)\\L 2 (n) yields the 
desired bound in the range s G [1, k + 2 — so]- □ 


4 FEM L 2 -error analysis on piecewise smooth geometries 

The convergence analysis of Theorem 13.31 did not make explicit use of the fact 
that a piecewise smooth geometry is considered; the essential ingredient was As¬ 
sumption 11.11 (which, of course, is related to the geometry of the problem). This 
is reflected in our use of 5r, which measures the distance from the boundary r. 
One interpretation of this procedure is that one assumes of the dual solution w 
(and, in fact, also of the solution of the “bidual” problem employed to estimate 
||<5 r 1 / 2+e (w — w/i)|| z, 2 (r 2 ) in Theorem l3.3l) that it may lose H 2 -regularity anywhere 
near r. However, for piecewise smooth geometries in conjunction with certain 
homogeneous boundary conditions (here: homogeneous Dirichlet conditions), this 
loss of H 2 -regularity is restricted to a much smaller set, namely, a subset of ver¬ 
tices in 2D and a subset of the skeleton (i.e., the union of vertices and edges) in 
3D. This set is given by M s in Definition 12.51 For this set M„, we introduce the 
distance function 

5 Ms :=dist(-,M s ), 5 M ,m=h-\-5 Mb . (4.1) 


Theorem 4.1 Let 17 be a polygon (in 2D) or a polyhedron (in 3D). Let M s be 
the set of vertices (in 2D) or edges and vertices (in 3D) associated with a loss of 
H 2 -regularity for 1 1. 3\) as given in Dehnition \2.5l Let Assumption ] 1.1\ be valid. Let 
Iu G Vh be arbitrary. Then we have 

||u - U h |U 2 (f2) < Ch\\6 S j^~ 1 'S/(u - Iu)\\ L 2 ( n y 
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Proof We may assume so < 1 since the case so = 1 corresponds to the standard 
duality argument with full elliptic regularity and set e := so — 1/2 G (0,1/2). The 
key observation is that, starting from the duality argument m , one can replace 
the weight function 5 r 1 ^ 2+e in (13.811 with any positive weight function. Taking as 
the weight function 5 jv/^ 2+£ , we get 

ll e lli 2 (j?) < II?m/ 2+£ V(u - In) |U 2 ( n ) ||4 / !“ £ V(« - w h )\\ LH ay (4.2) 

The estimate of w — Wh in the weighted norm is done similarly as in the proof of 
Theorem 13.31 taking into account the improved knowledge of the regularity of w. 
With SM s ,ch '■= {x G 17 | Sm 3 (x) < ch} we have the trivial bound 

II^m- s 2_£ V(w — u>h)\\L 2 (s 7 ) (4.3) 

^ \\^M 2 ~ e ^( W - W h)\\L^S Ms , oh ) + ll<^ s 2 “ £ V(w - W/ 1 )|U 2 (f2\S M „, eh )> 

where the parameter c will be selected sufficiently large below. The first term in 
(TOll is estimated in exactly the same way as in (13.1011 and produces 

¥m 2 ~ E V(w - Wh)|U 2 (s Msi0h ) < Ch\\e\\ L ^ n) . 

The second term in (TOll again requires the techniques from the local error anal¬ 
ysis of the FEM, this time with the appropriate modifications to account for the 
boundary conditions. Inspection of the arguments in |24l Sec. 5.3] shows that the 
key estimate (13.1111 extends up to the boundary in the following sense: 

||V(w - Wh)\\L*(B r c\n) 1$ l|V(w - Ihw)\\ L 2 ( Br , nr2 ) + — - ||w - w h \\L^( Br ,nn)\ 

(4.4) 

besides the implicit assumption r' > r + 0(h), the balls B r and B r ' are assumed 
to have the same center x and satisfy one of the following conditions: 

1. B r : = B r f(x) C 17; 

2. x G 917 and B r >(x ) n 17 is a half-disk; 

3. a: is a vertex of 17; 

4. (only for d = 3) x lies on an edge e and B r i (x) H 17 is a dihedral angle (i.e., the 
intersection of d(B r /(x) fl 17) with 917 is contained in the two faces that share 
the edge e. 

The reason for the restriction of the location of the centers of the balls is that 
the procedure presented in (241 Sec. 5.3] relies on Poincare inequalities so that 
the number of possible shapes for the intersections B r ' fl 17 should be finite. A 
covering argument (see Theorem IA.5I for the 2D case and Theorem IA.6I for the 
3D situation) then leads to the following bound with an appropriate ci > 0 (here, 
c > 0 is implicitly assumed sufficiently large): 

II^M 2_£V ( w — w h)\\L 2 (0\3 MatCh ) ^ (4.5) 

II^M S 2 ~ v ( w - I h w )\\L2(n\s M3 , cclh ) + II &M a /2 £ (w - w h )\\ L 2(n\s Me ,ec lh )- 
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The first term in (TOl) can be estimated with the improved regularity assertion of 
Lemma l2.6l to produce (with appropriate C 2 > 0 and the implicit assumption on c 
that CC 1 C 2 > 2) 

ll^f“ £ V( W - I k h w)\\ LH a\s M ^ lh ) < hfd^-^wW LHn\s Mmh) 

[ll^M s 2 ^lU 2 ^) + IMIff 3 / 2 +'(l?)] < h\\e\\ L 2( n y 

For the second term in (14.51) we note that —1/2 — e < 0 so that 5f^^ 2 ~ e < bf}^ 2 ~ e . 
This leads to 

II ^m] / 2 S ( wm w h)\\L2(n\s Ms ,oc lh ) ~\\Sm/ 2 £ ( w - w h )\\ L 2 (n) 

< ||<^ 1/2_e (w - w h )\\ L 2(n) < h~ 2E \\5~f 1/2+e (w - w h )\\ L 2 (n) ; 

the term h~ 2s \\5f, 1 ^ 2+e {w — Wh) ||l 2 (j 2 ) has already been estimated in (13.131) in the 
desired form. □ 

The regularity requirements on the solution u can still be slightly weakened. As 
written, the exponent so — 1 is related to the global regularity of the dual solution 
w. However, the developments above show that a local lack of full regularity of the 
dual solution w (and the bidual solution) needs to be offset by additional local 
regularity of the solution. To be more specific, we restrict our attention now to 
the 2D Laplacian, i.e., A = Id. In this case, the situation can be expressed as 
follows with the aid of the singular exponents ay := n/uij, where ay G (7t,27t) is 
the interior angle at the reentrant vertices Aj. j = 1, ..., J. 

Corollary 4.2 Let J7 C R 2 be a polygon and let A = Id. Let Sj := dist(-,Aj), 
j = 1 for the J reentrant corners. Set 5j := h + 5j. Let uy be the interior 

angle at Aj and aj = 7 r/wy. Fix /3y >1 — aj arbitrary. Then for any Iu G Vh 

J 

\\u-u h \\ L 2 {a) < l|£” ft V(u- Iu)\\ L 2 {n) . 

i =1 

Proof The proof follows by an inspection of how the regularity of the solution 
w = Te of the dual problem enters the proof of Theorem 14.11 By, e.g., (TO] the 
solution w = Te is in a weighted H 2 -space with 

\\Y[S?V 2 w\\L>M<\\e\\mn), ( 4 - 6 ) 

j =1 

and Assumption 11.11 holds with any so < minj ay. The regularity assertion (l4~6l) 
suggests to choose Yli= i ^ as the weight in the proof of Theorem [331 Inspection 
of the procedure in the proof of Theorem 14.II then leads to the result. □ 

We extract from this result another corollary that we will prove useful in the 
numerical results. We formulate it in terms of (standard, unweighted) Sobolev 
regularity in order to emphasize the difference in regularity requirements of the 
solution near the reentrant corners and away from them: 
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Corollary 4.3 Assume the hypotheses of Corollary \ f.2\ Let s > 1 and Si > 1, 
i — 1,..., J. Let U := 17 \ U Mi, for some neighborhoods Mi of the reentrant vertices 
Ai. Let u E H Si fUi), i = 1,..., J and u E H S {U). Then for arbitrary s > 0 

||u - uh\\l 2 {Q) < C e h T , r := min(l + fc, s, _min (-1 + aj + Sj — e)). 

Proof The approximant Iu in Corollary 14.21 may be taken as any standard nodal 
interpolant or the Scott-Zhang projection. Then standard estimates and Corol¬ 
lary 22] produce with the choice /3j := 1 — aj + e for arbitrary small but fixed 
£ > 0 : 

\\u-u h \\ L \n) <h min {/i mm{M_1} , h~ 0i+Si ~ 1 } 

< min 


□ 


5 Optimal T 2 (S'h)-convergence 

Additional regularity of the solution also allows us to prove that the error on the 
strip Sh of width 0{h) near T is of higher order: 

Theorem 5.1 Let Assumvtion \ 1. 1\ be valid. Then the FEM error u — Uh satisfies 
l|u - u h \\ L 2 ( Sh) < h k+3/2 ( 1 + <5fc,i| lnh\)\\u\\ B k+3/2 {nr 

\\u~ u h \\ L 2 {Sh) </i s+3/2 (l + 4 i i|ln/i|)||u|| B «+3/ 2(f2) , sE (0, k), 

||M - Uh\\Li{S h ) < /i 3/2 (l + 4,i|ln/i|)||u|| B 3/2 (r2) , 

where 5k,i is the Kronecker symbol. The implies constant depends on the shape 
regularity of the triangulation, 17, and the coefficient A. Specifically, it depends on 
ao and ||-A.11c°-i(T2) anc ^; the case k > 1, additionally on 11A11ci.3-( 72 ) - 

Remark 5.2 1. The regularity requirement (17) can be weakened: it suffices 

that u be in b\ +3 ^ 2 (Sd) in a fixed neighborhood Sd of T and in H k+1 (S 7). 
See [ 16| for the details of a closely related problem. 

2. Since (17) D i7 s+3 ^ 2 (l7), the assertions for s E (0 ,k) can be weakened 

by replacing ||u|| B »+3/2^^ with ||u||^f ii + 3 / 2 (r 2 ) on fi 16 right-hand side. Only for 
the limiting cases s = 0 and s = k, we require the stronger requirement u E 

B 2A 3 / 2 (tt) C H s+ 3 / 2 {n). . 

Proof The structure of the proof is very similar to that of Theorem 13.31 The 
main difference arises from the fact that the right-hand side of the dual problem 
is supported by the thin neighborhood Sh, and this support property has to be 
exploited. 
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Let e = u — Uh- Let xs h be the characteristic function of Sh- Let w = T(xs h e ) 
and Wh G Vh its Galerkin approximation. Again, Galerkin orthogonality for u — Uh 
and w — Wh implies 

||e||z, 2 (s h ) = (e,xe) = a(e,w) = a(e,w - w h ) = a(u - I%u,w - w h ) 

< C\\S~ 1/2 - £ V(u- ^u)\\L H o)\\5r /2+s V(w - w h )\\ LH a), ( 5 . 1 ) 


where e > 0 is arbitrary (in fact, e€l would be admissible). We flag at this point 
already that we will select e = 0 for k = 1 and e > 0 arbitrary (but sufficiently 
small) for k > 1. Each of the two factors in H5.ll) is estimated separately. 

1. step: For the first factor in G3D we use approximation properties of the 
Scott-Zhang operator together with Lemma l2.1l to get for j G {0,..., k} 


||«5- 1/2 - £ V(u 


ltu)\\ LHa) < V\\5- 1/2 - e V +1 u\\ L 2 ia) , 


f|ln/i| 1 / 2 ||V J ' + 1 u|| B i/a (fi) 

\h- £ \\V j+1 u\\ <hn) 


if e = 0 
if £ > 0. 


(5.2) 


With the Kronecker symbol <5o, £ , we have shown for j G {0,1,..., k} 

|| 5p 1/2 ~~\7{u - Ih.u)\\ L 2 (n) < h J h~ e ( 1 + 5 0 ,e| ln/i| 1/2 )||u|| B j+3/2 (f2) . (5.3) 


Since the Scott-Zhang operator is defined on 7L 1 (f2) irrespective of boundary 
conditions, we may use an interpolation argument to lift the restriction to integer 
values j. Specifically, the reiteration theorem (cf., e.g., |21[ Thm. 23.6]) asserts 
that the Besov space B^^^ 2 , which we have defined by interpolation between 

(integer order) Sobolev spaces, coincides with the interpolation space between 
Besov spaces, viz., 

B 2 + 1 /2 W) = ( B l,[i (12), B 2 x 3/2 (ft)) s /k,oo (equivalent norms). 


Hence, we may decompose for arbitrary t > 0 a function u G B%\£J 2 (Q), s G (0, fc), 

as u = u — u\ + ui with ui G B^ 3 ^ 2 (17) and uo '■= u — u\ G B^ 2 (fl) together 
with 


Ml B 3 2 '?(n) ^ Ct s/k \\u\\ 


H Wi Hb2+ 3/2 (17) — Ct 


i/fc-1. 


s+3/2 


W 


Writing u — I^u = (uo — Ih u o) + (ui — Z^ui) we can use (15.311 with j = k for the 
second term in brackets and j = 0 for the first term in brackets to get with the 
choice t = h k 


\\5 r 1/2 c V(u — Ihu)\\ L 2 (n) < h s h e (l + 5o,e| ln/i| 1/2 )||u|| B »+3/2 (fi) . (5.4) 

Combining the estimates (15.31) with j = 0 and j = k and (15.41) for s G (0, k) we 
arrive at 


||^ 1/2 - 8 V(u-liu)|| L 2 (fl) 


< (1 + 5 0 , e | In h| 1/2 )/i _£ 


^ H u IIb“' i ; 3 / 2 (i?)’ 

h II m IIbJ+ 3 / 2 (i7)’ 


s = 0, 
s G (0, fc), 
s = k. 


(5.5) 
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2. step: The second factor in m requires more work. We start with a regu¬ 
larity assertion for w that exploits the support properties of xs h e an d follows from 
(ITTTD and (15T2)l : 

lhll B 3/^ (fi ) )$ h 1/2 \\xs h e\\mn), (5.6) 

IMIff3/2 +e(fi) < h 1/2 ~ £ \\xs H e\\ L 2 (n), e £ (0, s 0 - 1/2]. (5.7) 

We obtain an energy error estimate for w — Wh in the standard way by using 
quasi-optimality, the approximation properties of Vh, and the regularity assertion 

(EDI : 

Ik - w h \\ H i (n) < Ik - £ ^ 1/2 |kll B ^ ( fi) ~ h \\xs h e\\ L 2 ( o). (5.8) 

Lemma [3.1 1 is applicable with z = w — Wh', hence, obtain with (15.811 

W 5 r 1/2+e (™ ~ w h )\\mn) < h 3/2+e \\xs h e\\ L 2 (n), e G (0, s 0 - 1/2], (5.9) 
\\5r 1/2 (w - w h )\\ L 2(n) < h 3/2 \ ln/i| 1/2 ||xs Jl e|| Il 2 (fi) . (5.10) 

The bound (ED) informs us that control of w — Wh in a weighted H 1 -norm is re¬ 
quired. In this direction, we first write for a constant c > 0 that will be determined 
later sufficiently large 

\\5 1 I ! 2+e V (w — Wh)\\L 2 (S7) 

< \\5 1 I ! 2+e V(w — 'iiift,)||i2 ( s oh) + ||^ /2+£ V(w- w h )\\ L ^ n \ Sch ) 

< Ch 1/2+£ \\S7(w - w h )\\ LHn) + ||^ /2+£ V( W - w h )IU* ( o\s.H ) 

? Ch 3 / 2+e || XSft e|U2 (fi) + ||^ 2+£ V(u;-^)|U 2(fi \ Sch) . (5.11) 

We emphasize that e = 0 is allowed in (15. 1111 . It remains to control ||(5^/ 2+c V(u) — 
Wh)\\L 2 (n\s ah )- This is done again with the same arguments from the local error 
analysis as in the proof of Theorem ED The estimate (13.1211 holds verbatim, i.e., 

ll4 /2+e V( W - Wh)\\ L *(n\s ch ) (5.12) 

< ||4 /2+S V( W - lHw)\\ LHO \s cclh ) + ll^ 1/2+£ k - Wfc)IUw« lfc )- 

We emphasize that e = 0 is admissible in (13.1211 . As in the proof of Theorem 13.31 
the constant c will be selected in dependence of various inverse estimates that are 
applied. Combining (15.911 . (15.101) . (15.111) . (15.1211 we see that we have shown 

\\5 1 r /2+£ S7(w - w h )\\L*(n) (5.13) 

< ih 3/2+e \\xs h e\\mn) + ||^ /2+£ V(u; if £ > °> 

~ k 3/2 | ln/i| 1/2 ||xs,e|U2 (f2) + ||^ /2 V( W - Jjt«;)|| i 2 ( „ xgcich) if e = 0. 

3. step: We estimate the approximation error ||<5]k +c V(u> — 7h w )llz. 2 (f3\S 0 ich)- At 
this point the cases k = 1 and k > 1 diverge: since w solves a homogeneous elliptic 
equation on J7 \ S cic h (if cic > 1), interior regularity is available so that higher 
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order approximation can be brought to bear if k > 1 in contrast to the case k = 1 . 
We start with the simpler case k = 1. 

The case k = 1: From standard approximation results for I k , the inverse esti¬ 
mate of Lemma l2.9l (0, and (15.611 we get for a constant c 2 G (0,1) (implicitly, we 
assume that c is so large that c 2 cic/i > 2 h) 

ll^ / 2 V( W - l£w)\\ LH n\ Scich ) < 

< h\\nh\ 1/2 \\w\\ B 3/^ {n] < /i 3 / 2 |ln/i| 1 / 2 ||xs h e|| I , 2(fi) . (5.14) 

Inserting (15.51) (with e = 0) with the combination of (15.1411 and (15.1311 (again 
with e = 0) in (15.111 yields the desired final estimate for the case fc = 1 if we fix 
c = 2/(ciC 2 ). 

The case k > 1; We fix an e G (0,so — 1/2] arbitrary. From standard ap¬ 
proximation results for I k , the inverse estimates of Lemma [2.91 and the regularity 
assertion «£Z> we get (again for suitable constants c 2 , C 3 G (0,1) and the implicit 
assumption that c is such that C 3 C 2 cic is sufficiently large) 


||^ 2+£ V( W -7^)|| I , 2( ^ Scich) < h 2 ||^ /2+£ V 3 HU 2 (rA*WO 

Lem. [2T9] (fzTl r _ . 

< h 2 \\\6r 1/2+ ^ 2 W \\L H n\s mch) + +£ Vw\\ LWmoh ) 


< h 2 1/2+E [l|V 2 w|| L 2 (17 \ Sc3C2cioh) + ||Vw|| L 2 (r2 \ Sc3C2cich) j 

Lem. 12.91 \ivl 


< 


h 


2-1/2+ei 


Xs h e\\ L 2(n)- 


(5.15) 


Combining this with (15.511 produces in (15.111 the desired final estimate for the case 
k > 1 . □ 


From Theorem 1 5.1 1 we can extract optimal convergence estimates for the flux error 
I! dn(u - Uh)\\L*(ry- 

Corollary 5.3 Let Assumvtion \l. l\ be valid. Then with the Kronecker symbol 5 j-i 


II d n u - d n u h \\ L 2 (r) < (1 + 4 ,i| lnh|) < 


h ll' u ll B fc+3/ 2 ( f 2)! 

^ II^11 ^ ^ (0?^): 

. II^H ' 


Proof Structurally, the proof follows jl8l Cor. 6.1] in that estimating the error on 
T is transferred to an estimate on the strip Sh- The triangle inequality gives 

\\dn(u-u h )\\ LHn < II dn{u - Ihu)\\ L 2 ( r) + || d„(IhU - u h ) IIi 2 ( r). (5.16) 

The two terms in (15.1611 are estimated separately. 

1. step: We claim that 


{ h k \\u\ 


II dn(u - 7^u)|| L 2( r ) < < 


B k + 3 , \n) 


^ II^11 j ^ ^ 


(5.17) 
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We will only show the limiting cases u G B^ 3 ^ 2 (17) and u G B'^ 2 (£2)', the inter¬ 
mediate cases follow by an interpolation argument similar to the one used in the 
proof of Theorem 15.11 For the case of maximal regularity, we use an elementwise 
multiplicative trace inequality for the elements abutting T to get 

II dn(u - I k h u )IMn < ft fc/ Vllv fc+ 1 «IU» ( s ah) ft (fc " 1 )/ Vnv fc+ 1 “ll^ ( ^ ) 

, J231 

<h k - 1/2 \\S7 k+1 u\\ LHS2h ) < /i fc ||V fe+1 M|| B i/ 2 (fi) < h k \\u\\ B k+3/i^ n y 

For the case of minimal regularity, u G B 3 J 2 ( fl) we first note that we obtain from 
(12.51) that ||u|| L 2 (jq < \\v\\ B i/ 2 ^ n y Using this and inverse estimates, we get 

\\d n (u — Ih u )\\L 2 {r) < \\d n u\\L2(r) + \\dnIh u \\L 2 (r) 

<\\Vu\\ Bl ^ n) + h ~ 1 / 2 \\VlU\LHS h ) 

/^stable 1231 

£ II^Hb^w + h " 7 llv«IU* ( s aJ .) < llv«|| B i/ a(n) . 

2. step: The term \\d n (IhU — Uh)\\L 2 (r) in (15.161) is controlled with inverse 
estimates and Theorem 15. II as follows: 

\\dn(Ih.u - u h )\\ L 2 (r) < h -1 2 \\V{Ihu - u h ) ||z,2(s h ) ^ h~ 3/2 \\IhU - u h \\ L 2 (s h ) 

< h~ 3 2 \\u - Ihu\\ L 2 {Sh ) + h~ 3/2 \\u - u h \\ L 2 (Sh) . 

The term \\u — I^u\\l 2 {s h ) can b e controlled with the approximation properties of 
Ih in the desired fashion: ||u - l£u\\ L 2 (Sh) < /i||Vu||ia (S2h) < h 3/2 ||Vu|| B i/a (fi) . 
The contribution ||w — Uh\\L 2 (S h ) is estimated with the aid of Theorem 15.II □ 


6 Extension of the results of [18| 

The arguments of the present paper are similar to those underlying [TH] , in spite of 
the fact that we did not employ the anisotropic norms that we introduced in |18| 
but instead worked with weighted Sobolev spaces. A feature of the analysis here 
that was not present in [18] is our FEM error analysis in Lemma r3.1l for a weighted 
L 2 -estimate, which, in turn, relies on the regularity assertions of Lemma 12.81 for 
problems with data in weighted spaces. This additional technical issue was cir¬ 
cumvented in 1 18] by assuming convexity of Q so that optimal order L 2 -estimates 
could be cited from the literature. The present analysis provides the necessary 
technical tools to remove this simplification in |18| , where a more complex mortar 
setting is analyzed. It is possible to make use of weighted L 2 -estimates similar to 
those of Lemma 13. II in the setting of |T8] , For that, regularity results of the type 
provided in Lemma f2.8l have to be used. The outcome of this refinement is that the 
main results of [18] . namely, [18[ Thm. 2.1], which provides L 2 -estimates on strips 
of width 0(h) around the skeleton, and [18l , Thm. 2.5], which provides optimal 
order approximations for the mortar variable, hold true if the geometry is such 
that Assumption ll.ll is valid. We will not provide the details of the arguments here 
and refer to El Appendix B] instead. Nevertheless, for future reference we record 
the end result: 
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Theorem 6.1 In [T8[ Thms. 2.1, 2.5], the assumption of convexity of ft can be 
replaced with \18 1 Assumption (5.2)]. 


7 Numerical results 

We consider the simple model equation — Au = f in ft C R, d G {2, 3} with inho¬ 
mogeneous Dirichlet boundary conditions. These are realized numerically by nodal 
interpolation of the prescribed exact solution u , and the data / is also computed 
from u. In the case of a non-smooth solution, we use a suitable quadrature formula 
on finer meshes to guarantee that the L 2 -error is accurately evaluated. 


7.1 Two-dimensional results 

We use a sequence of uniformly refined triangular meshes, where each element is 
split into four triangles. 


7.1.1 Lowest order discretization 

We consider two typical domains for reentrant corners. We start with the L- 
shaped domain (—1, l) 2 \ [0,1] x [—1,0] and then consider a slit domain (— 1, l) 2 \ 
((0,1) x {0}). In both cases, the prescribed solution is given in polar coordinates 
by u(r,(j>) = r“ sin(a/>) where a, a are given parameters. For non-integer a, we 
have u G B\ + ^(ft) by j4j Thm. 2.1]. Moreover, we test the influence of the position 
(xo,yo) of the weak singularity at r = 0 by defining r 2 := (x — xo) 2 + (y — Vo) 2 - We 
note that irrespective of the location ( xo,yo ) of the singularity on the boundary 
r, we have u G (ft) C H 1+a ~ e (ft) for any e > 0. 



(so, Vo) = (0,0) 

a = 7t/2 

(xo,yo) = (0.5,0) 

a = 7T 

(x 0 ,yo) = (0,1) 
a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

81 

6.1585e-03 

- 

6.8141e-03 

- 

6.2506e-03 

- 

289 

2.6986e-03 

1.19 

2.5648e-03 

1.41 

2.1211e-03 

1.56 

1.089 

1.1123e-03 

1.28 

8.8428e-04 

1.54 

6.7413e-04 

1.65 

4.225 

4.4037e-04 

1.34 

2.9202e-04 

1.60 

2.0903e-04 

1.69 

16.641 

1.7107e-04 

1.36 

9.4164e-05 

1.63 

6.4027e-05 

1.71 

66.049 

6.5689e-05 

1.38 

2.9909e-05 

1.65 

1.9471e-05 

1.72 

263.169 

2.5030e-05 

1.39 

9.4012e-06 

1.67 

5.8930e-06 

1.72 

1.050.625 

9.4877e-06 

1.40 

2.9328e-06 

1.68 

1.7774e-06 

1.73 

4.198.401 

3.5834e-06 

1.40 

9.0968e-07 

1.69 

5.3475e-07 

1.73 


Table 7.1 L-shaped domain, k = 1: Influence of the position of singularity for a = 0.75. 


For the L-shaped domain, the shift parameter so can be taken to be any so < 
2/3. From the theoretical results in Section [3j we therefore expect the error decay 
to have a rate of at least min(2,1+a —1/3) uniformly in the position ( xo , yo) of the 
singularity. Table [7~i1 shows the numerical results for a = 0.75 and a = 2/37T, for 
which min(2,1+a—1/3) = 1.417. As it can be seen for ( xo , yo) = (0, 0), we observe 
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a good agreement with Theorem l3.3l However for the locations (xo,yo) = (0.5,0) 
and (xo, yo) = (0,1), the rates are substantially better. This can be explained by 
the more refined analysis of Section}!] Using Corollarv l4.31 we expect an improved 
convergence rate of 1.75 for these cases. 



a = 10/9 

a = 4/3 

P 

II 

CO 

to 

DOFs 

2-error 

rate 

L 2 -error 

rate 

L z -error 

rate 

81 

6.5660e-03 

- 

8.6776e-03 

- 

8.9932e-03 

- 

289 

2.3309e-03 

1.49 

2.8523e-03 

1.61 

2.8151e-03 

1.68 

1.089 

7.3413e-04 

1.67 

8.2870e-04 

1.78 

7.8034e-04 

1.85 

4.225 

2.2257e-04 

1.72 

2.3073e-04 

1.84 

2.0751e-04 

1.91 

16.641 

6.5650e-05 

1.76 

6.2539e-05 

1.88 

5.3910e-05 

1.94 

66.049 

1.9056e-05 

1.78 

1.6688e-05 

1.91 

1.3835e-05 

1.96 

263.169 

5.4810e-06 

1.80 

4.4099e-06 

1.92 

3.5256e-06 

1.97 

1.050.625 

1.5690e-06 

1.80 

1.1580e-06 

1.93 

8.9467e-07 

1.98 

4.198.401 

4.4822e-07 

1.81 

3.0279e-07 

1.94 

2.2641e-07 

1.98 


Table 7.2 L-shaped domain, k = 1: Influence of exponent a for a = 2/3 n and (xo , ijq) = (0, 0). 


Table f7~2l shows the results for (xo,yo) = (0,0) and a E {10/9,4/3,3/2}. From 
Theorem 13.31 we expect convergence rates of 1.78, 2, and 2, respectively. The 
observed numerical rates of 1.81, 1.94, and 1.98 are quite close. 

The situation is similar for the slit domain where the regularity of the dual 
problem is even further reduced. It corresponds to a limiting case of our theory, 
which, strictly speaking, we did not cover, since the parameter so of Assumption ! 1. ll 
may be taken to be any so < 1/2. Nevertheless, one expects from Theorem 13.31 a 
convergence rate close to min{2,1 + a — 1/2}. For a = 0.75 this is 1.25, which is 
visible in Table rf~3i for the case (xo,yo) = (0,0). Again, the better convergence 
behavior for (xo,yo) = (0.5,0) and (xo,yo) = (0,1) can be explained by the the¬ 
ory of Corollary 231 which predicts 1 + a = 1.75. Table I7T1 shows the results for 
(xo,yo) = (0,0) and a E {10/9,4/3,3/2}. From Theorem 13.31 we expect conver¬ 
gence rates of 1.61, 1.83 and 2, respectively. The observed numerical rates of 1.65, 
1.86, and 1.96 are reasonably close to these predictions. 



Oo,yo) = (o,o) 

a = 7r/2 

{x 0 ,yo) = (0.5,0) 
a — 7r 

{x 0 ,yo) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

6.1391e-03 

- 

1.1088e-02 

- 

1.0692e-02 

- 

348 

2.8187e-03 

1.12 

4.1329e-03 

1.42 

3.8553e-03 

1.47 

1.315 

1.2351e-03 

1.19 

1.4164e-03 

1.54 

1.3388e-03 

1.53 

5.109 

5.3338e-04 

1.21 

4.7830e-04 

1.57 

4.4562e-04 

1.59 

20.137 

2.2846e-04 

1.22 

1.4725e-04 

1.70 

1.4420e-04 

1.63 

79.953 

9.7267e-05 

1.23 

4.6683e-05 

1.66 

4.5843e-05 

1.65 

318.625 

4.1233e-05 

1.24 

1.4761e-05 

1.66 

1.4401e-05 

1.67 

1.272.129 

1.7428e-05 

1.24 

4.3773e-06 

1.75 

4.4861e-06 

1.68 

5.083.777 

7.3524e-06 

1.25 

1.3285e-06 

1.72 

1.3889e-06 

1.69 


Table 7.3 Slit domain, k = 1: Influence of the position of singularity for a = 0.75. 
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a = 10/9 

a = 4/3 

a = 3/2 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

5.7534e-03 

- 

7.3549e-03 

- 

7.5901e-03 

- 

348 

1.9412e-03 

1.57 

2.2414e-03 

1.71 

2.1664e-03 

1.81 

1.315 

6.2583e-04 

1.63 

6.4849e-04 

1.79 

5.8638e-04 

1.89 

5.109 

1.9689e-04 

1.67 

1.8251e-04 

1.83 

1.5450e-04 

1.92 

20.137 

6.1446e-05 

1.68 

5.0718e-05 

1.85 

4.0197e-05 

1.94 

79.953 

1.9191e-05 

1.68 

1.4021e-05 

1.85 

1.0396e-05 

1.95 

318.625 

6.0229e-06 

1.67 

3.8699e-06 

1.86 

2.6803e-06 

1.96 

1.272.129 

1.9023e-06 

1.66 

1.0682e-06 

1.86 

6.8978e-07 

1.96 

5.083.777 

6.0474e-07 

1.65 

2.9514e-07 

1.86 

1.7730e-07 

1.96 


Table 7.4 Slit domain, k= 1: Influence of exponent a for a = 1 /2tt and (xq, yo) = (0,0) 


7.1.2 Second, order finite elements 

In this subsection, we test the performance of quadratic finite elements for the L- 
shaped domain. We use the same type of solution as before and vary the parameter 
a for (xo, yo) = (0, 0), i.e., the re-entrant corner. Here we expect from our theory a 
convergence rate of min(3, a + 1 — 1/3). For a G {2.175,2.275, 2.375}, the observed 
numerical rates, which are visible in Table 1731 are very close to the theoretically 
predicted ones. 



a = 2.175 

a = 2.275 

a = 2.375 

DOFs 

L2 error 

rate 

L2 error 

rate 

L2 error 

rate 

289 

2.7565e-04 

- 

2.4570e-04 

- 

2.2177e-04 

- 

1.089 

5.1121e-05 

2.43 

4.1696e-05 

2.56 

3.3912e-05 

2.71 

4.225 

7.5320e-06 

2.76 

5.7319e-06 

2.86 

4.3221e-06 

2.97 

16.641 

1.1051e-06 

2.77 

7.8407e-07 

2.87 

5.4888e-07 

2.98 

66.049 

1.5938e-07 

2.79 

1.0553e-07 

2.89 

6.8762e-08 

3.00 

263.169 

2.2723e-08 

2.81 

1.4044e-08 

2.91 

8.5292e-09 

3.01 

1.050.625 

3.2138e-09 

2.82 

1.8538e-09 

2.92 

1.0497e-09 

3.02 


Table 7.5 L-shaped domain, k = 2: Influence of a for a = 2/3n and (xq , yo) = (0,0). 


7.2 Three-dimensional results 

In the three dimensional setting, we consider a Fichera comer fl := (—1,1) 3 \ [0, l] 3 
and prescribe the smooth solution u(x,y,z) := sin((x + y)7f) cos(2nz). The inho¬ 
mogeneous Dirichlet conditions are realized by nodal interpolation. The discretiza¬ 
tion is based on trilinear finite elements on hexahedra and uniform refinements. 
Although the dual problem lacks full regularity, Theorem 13.31 asserts that this 
can be compensated by extra so regularity of the primal solution to maintain full 
second order convergence in L 2 . 

Table ITTfil shows that we observe numerically already for coarse discretizations 
the predicted convergence order two, and the theoretical results are confirmed. 

We point the reader to Appendix [E] for a further numerical results. 
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DOF 

L z -error 

rate 

316 

3.032 

26.416 

220.256 

1.798.336 

14.532.992 

0.075444 

0.017182 

0.0039376 

0.00094597 

0.00023208 

5.7491e-05 

1.96 

2.04 

2.02 

2.01 

2.00 


Table 7.6 Fichera corner, k = 1: L 2 -error for a smooth solution. 
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A Coverings 

In this appendix, the distance dist(a:,M) for some set M appears frequently. For 
notational convenience, we set dist(a;, 0) = 1 to include the degenerate case M = 0. 
We quote from |18l Lemma A.l]: 

Lemma A.l Let 12 C R d be bounded open and M = M be a closed set. Fix 
c G (0, 1) and £ G (0, 1) such that 

l-c(l + e) =:c 0 >0. (A.l) 

For each x G 12, let B x := B c dist(x,M)( x ) be the closed ball of radius cdist(x,M) 
centered at x, and let B x := -Bjq+e^dist^.M)^) denote the stretched (closed) ball 
of radius (1 + e)cdist(a:, M ) also centered at x. 

Then there exists a countable set Xi G 12, i G N, and a constant N G N 
depending solely on the spatial dimension d with the following properties: 

1. (covering property) Ui^B Xi D 12; 

2. (finite overlap on f2) for each x G 12, there holds card{i | x G B Xi } < N. 

Proof [181 Lemma A.l] assumed that M C 12. However, an inspection of the proof 
shows that this is not necessary. 

Before we proceed with variants of the covering result of Lemma rA.il we introduce 
the notation of sectorial neighborhoods relative a singular set M: 

Definition A.2 (sectorial neighborhood) Let e, M Cl d and c > 0. Then 

S e ,M,c := Ux6e23Hdist(a;,M) {x) 

is a sectorial neighborhood of the set e relative to the singular set M. 

We are interested in coverings of lower-dimensional manifolds by balls whose cen¬ 
ters are located on these manifolds: 

Lemma A. 3 Let d G N and 1 < d' < d. Let u> C R d and let 12 C R d be the 
canonical embedding of to into R d , i.e., 12 := lj x {0} x • • • x {0} C R d . Assume 
the hypotheses and notation of Lemma \A.l\ Then there are c > 0, N > 0, and a 
collection of balls B Xi , iGN, as described in Lemma | A. 1 1 such that 

(i) (covering property for L2) Ui^B Xi D 12. 

(ii) (covering property for a sectorial neighborhood of L2) U i£®B Xi D Sq^m,c- 
(Hi) (finite overlap property on for each x G R d , there holds card{i | x G 
B Xi } < N. 

Proof We employ the result of Lemma IA.1I for the lower-dimensional manifold w 
noting that B x fi u> is a ball in R d . In order to be able to ensure the covering 
condition for the sectorial neighborhood of 12 stated in lHull . we introduce the 
auxiliary balls B' x := -B c /2dist(x,M)(20 of half the radius. Applying Lemma lA.ll 
with these balls B x and the stretched balls B x therefore produces a collection of 
centers Xi G 12, i G N, such that 

1. B' Xi n 12 covers 12; 
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2. for the stretched balls B Xi , we have a finite overlap property on 12: 

Vacard{i | x G B Xi } < Al. (A.2) 

We next see that the balls B Xi even have the following, stronger finite overlap 
property: 

Va; G R d : card{i | x G B Xi } < N. (A.3) 

To see this, define the infinite cylinders C Xi := {x \ TTd'(x) G B Xi n 12}, where ix# 
is the canonical projection onto the hyperplane {x = (xi,..., Xd) G R d | Xd’+i = 

■ ■ ■ = Xd = 0}. Clearly, B Xi C C Xi . These infinite cylinders have a finite overlap 
property by (TOll as can be seen by writing any x G R d in the form x = (nd' ( x ), x') 
for some x’ G R d-d and then noting that x G C Xi implies ■Kd'ix) G B Xi fl 12. 

Is remains to see that the balls B Xi cover a sectorial neighborhood of 12. To 
that end, we note that the balls B x . cover 12. Furthermore, for each x G 12, 
we pick Xi such that x G B x . C B Xi . Since the radius of B Xi is twice that of 
B x ., we even have B c / 2 dist{x i ,M){ x ) C B Xi . Furthermore, by c G (0,1), we have 
0 < (1 — c/2) dist(xi, M) < dist(x, M) < (1 + c/2) dist(xi, M). Therefore, there is 
c > 0 such that B^dist(x,M)( x ) C B Xi and thus 

U x ^qBc dist(rc,M) ( x ) c U iB Xi . 


□ 

We next show covering theorems for polygons and polyhedra. In the interest of 
clarity of presentation, we formulate two separate results. Before doing so, we point 
out that balls with center located on the boundary of the polygon/polyhedron 12 
will feature importantly so that the intersection of this ball with 12 will be of 
interest. We therefore introduce the following notions: 

Definition A. 4 (solid angles and dihedral angles) 

1. Let 12 C R 2 be a Lipschitz polygon. Let A be a vertex where the edges ei, ei 
meet. We say that the set B e {A) n 12 is a solid angle, if d(B e (A) fl 12) fl 912 is 
contained in {A} U ei U e 2 - 

2. Let 12 C R 3 be a Lipschitz polyhedron. Let A be a vertex of 12. We say that 
the set B e (A) (~l 12 is a solid angle , if d(B e (A ) fl 12) fl 912 is contained in the 
union of {A} and the edges and faces meeting at A. 

3. Let 12 C R 3 be a Lipschitz polyhedron. Let e be an edge of 12, which is shared 
by the faces /i, /y. Let x G e. We say that the set B e (x) fl 12 is a dihedral angle, 
if d(B e (x) fl 12) fl 912 is contained in eU/iU/ 2 . 

Theorem A.5 Let 12 C R 2 be a bounded Lipschitz polygon with vertices Aj, j = 
1,..., J, and edges £. Let M C {Ai,..., Aj}. Set A' := {Ai,..., Aj} \ M and fix 
£G (0,1). 

(i) There is a sectorial neighborhood °f the vertices A' and a constant 

c G (0,1) such that Sjp^m,c covered by balls Bi := B c dist( Xi ,M)( x i ) with 
centers Xi G A' . Furthermore, the stretched balls Bi := fi(i+ e ) c dist(rc i ,M)(2 ; i) 
are solid angles and satisfy a finite overlap property on R 2 . 
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(ii) Fix a sectorial neighborhood IA := Sa',m,c' of the vertices A'. For each edge 
e G £, there is a sectorial neighborhood S ei M,c an d a constant c G (0,1) such 
that S e ,M,c\M is covered by balls Bi = -B c dist(xi,M) {x%) whose centers Xi 
are located on e. Furthermore, the stretched balls Bi = B( i+ e ) c dist(xi,M)(*i) 
satisfy a finite overlap property on R 2 and are such that each Bi Pi 17 is a 
half-disk. 

(in) Fix a sectorial neighboodlA := > Ss t M,c' of the edges £. There is c G (0,1) such 
that £2\IA is covered by balls Bi = -Bcdistlxi.M)^*) such that the stretched 
balls Bi = -B(i_|_ e ) c dist(xi,M) (xi) are completely contained in £2 and satisfy a 
finite overlap property on R 2 . 

Proof The assertion (0 is almost trivial and only included to emphasize the struc¬ 
ture of the arguments. Assertions ©, (1ml) follow from suitable applications of 
Lemmas IA.3I and IA.1I □ 

The 3D variant of Theorem lA.5l is formulated in Theorem IA.61 We emphasize that 
the “singular” set M need not be the union of all edges and vertices but can be 
just a subset. We also emphasize that it is not necessarily related to the notion of 
“singular set” in Definition l2.5l although it is used in this way. The key property of 
the covering balls is again such that the centers are either a) in 17 (in which case 
the stretched ball is contained in 17); or b) on a face (in which case the stretched 
ball Bi is such that Bi fl 17 is a half-ball); or c) on an edge in which case Bi D 17 
is a dihedral angle (see Definition IA. 41) : or d) in a vertex in which case Bi fl 17 is 
a solid angle (see Definition IA.41) . 

Theorem A.6 Let 17 C R 3 be a Lipschitz polyhedron with faces T, edges £, and 
vertices A. Let M C A and Mg C £. Let M = M = M _4 U Mg and fix e G (0,1). 
Let A' := {A G A \ A £ M} be the vertices not in M and £' := (e G £ | efl M = 0} 
be the edges not abutting M. Then: 

(i) (non-singular vertices) There is a sectorial neighborhood Sa',m,c of the ver¬ 
tices in A! and a constant c G (0,1) such that S^ t M,c is covered by balls 
Bi := B c dist(a; i ,M) ( x i) with centers Xi G A! . Furthermore, the stretched balls 
Bi := B( i+ e )cdist(xi,M) (**) are solid angles and satisfy a finite overlap prop¬ 
erty on R 3 . 

(ii) (non-singular edges) Fix a sectorial neighborhood Li := Sa',m,c' of A'. For 
each edge e G £', there is a sectorial neighborhood S e ^m,z an d a constant 
c G (0,1) such that S Et M,c \ M is covered by balls Bi = B c dist)^, M){xi) 
whose centers Xi are located on e. Furthermore, the stretched balls Bi = 
-B(i_l_ e ) c dist(:Ej,M) (*£i) satisfy a finite overlap property on R 3 and Bi D 17 is a 
dihedral angle. 

(Hi) (faces) Fix a sectorial neighboodlA := Se.,M,d of £. There is a sectorial neigh¬ 
borhood Sjr,M,c and a constant c G (0,1) such that Sjr tM ^\U is covered by 
balls Bi = B c dist(xi ,with centers Xi G 917. Furthermore, the stretched 
balls Bi = S(i+ e )cdist(a: i) M)(^i) satisfy a finite overlap property on R 3 and 
Bi Pi 17 is a half-ball. 

(iv) (interior) Fix a sectorial neighbood IA := of T, where T is the set 

of faces. Then there is c G (0,1) such that £2\U is covered by balls Bi = 
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B c dist(xi,M)( x i) with centers Xi G 1?. Furthermore, the stretched balls Bi = 
-®(i+e)cdist(a:i,M) ( x i) satisfy a finite overlap property on R 3 and Bi C fl. 

Proof Follows from Lemmas IA.3I and IA.1I □ 
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B Details for the extension of the results of [T8j . 

The arguments used in the proof of Theorem 15.11 rely on techniques developed in 
jl~8] . A feature of the analysis here that was not present in m is the bidual problem 
with right-hand side 5f 1 (w — Wh ) that allowed us to estimate w — Wh in a weighted 
space. This technical issue was circumvented in ||18| by assuming convexity of 
fi so that optimal order L-estimates could be cited from the literature. The 
corresponding bidual problem can be analyzed in the mortar setting of [18] as well. 
The end result is then Theorem iBTl which states that the convexity assumption 
in [18] can be relaxed to the validity of Assumption 11.11 for the Poisson problem, 
i.e., [HI (5.2)]. 

In the interest of brevity, we employ in this appendix the notation of |18] and 
assume the reader’s familiarity with ]18| . 

It will be useful to write Hf iw (r) for the space given by the broken Sobolev 
norm on the skeleton T, i.e., the Sobolev norm is understood facewise. Further¬ 
more, we will write || ■ || ff i for the broken ff 1 -norm, i.e., || ■ ||^i = JT || ■ IIh 1 ^)- We 

also introduce the L 2 -projection 77^, '■ L 2 {r) —> Mh and recall that, since Mh is 
a product space based on the faces, it inherits from j!8l (A2)] the approximation 
property 

\\z - nM h z\\ L 2 (r) < Ch s \\z\\ H s^ r) , s G [0, k\. (B.l) 

The main result is: 

Theorem B.l In Jl8i Thms. 2.1, 2.5], the assumption of convexity of Q can be 
replaced with HE Assumption (5.2)]. 

Proof We will only sketch the modifications entailed by the weakenend regularity 
assumptions. 

Proof of [18l Thm. 2.1]: The starting point is the error representation [181 
(6.2)], which consists of three terms: 

a(w - w h ,u- P h u) + b(w - w h , A - p h ) + b(u - P h u,\ w - p h ), (B.2) 

where p-h, ph £ Mh are arbitrary. The first term in (IB.211 can be estimated as in [l8l, 
Proof of Thm. 2.1] in view of the generalization of [T8I Lemma 5.5] given below as 
Lemma [B.3l The third term in (IB.21) can again be estimated as in [181 (6.4)] since 
Lemma |B.2l below provides the estimate inf,-. <= w. ||Am — Mh.||.L 2 (r) ^ 

The second term in (lB~2ll requires a modification of the procedure in [ 18l (6.3)]. 
Taking ph = n]f h X in [18] (6.3)] yields 

b(w - w h , A - IlM h \) = I [w - w h ](\ - IlM h \) 

= I ([to - Wh] - nM h [w - iufc])(A - J7mia) 

< h 1 '- 1| [w - w/,]|| ff i/2 (r) ||A - I7M h ||i 2 (r) 

< h 1 ' ~\\w — Wh\\m IIA — nM h M\L 2 (r) 

< h 1/2 1 k \\w - Wk||ffi||A|| 
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the last step followed from (IB. Ill . The proof is completed with the aid of Lemma[h2] 
and the trace estimate ||A||#k (jq < \\u\\ B k+3/2^y 

Proof of |18l Thm. 2.5]: The proof stands as given in |18] . □ 

Lemma B.2 Assume that Q satisfies J181 Assumption (5.2)]. Then, for v G 
L 2 (Sh) C L 2 (17) and w := T D (v) and the corresponding Lagrange multiplier X w 
defined facewise by Am| 71 = —9 n w|f2 s(i) and the corresponding mortar approxima¬ 
tion u>h of w there holds: 



w h || 


2 H i( fii ) = W w ~ w h\\m < h\\v\\ L 2 {n) , (B.3) 

||A™ — L[M h \ w \\ L 2( r ' ) < /i 1 ' /2 ||u||L 2 (r2)- (B.4) 


Proof We start with the proof of (IB. 31) . It results from standard convergence theory 
for mortar methods as follows. [181 Assumption (5.2)] provides w G H' i ^ 2+£ [T2) for 
some £ > 0 together with ||ui|| H 3/2+ c ^ < C |MIh-i/ 2 +'(^)' The L a g ran g e multi¬ 
plier Au, is given facewise by the expression Au,| 7i = —d n w\n sm G H e (yi) together 
with the estimate ||A™||#e (r) ||n||The standard convergence theory 

for mortar methods (as worked out, e.g., in [8] Prop. 2.3]) then gives 


Ik “ w h \\ H i < h 1/2+1 \\w\\ H 3,2 +c{n) < h 1/2+£ \\v\\ H - 1/2+c{n) . 
The proof of (|B.3ll is complete if we can show that 


IMIff-i/2+c(fi) < h 1/2 : |M| L 2 (ifi) . (B.5) 

This last estimate exploits supp v C Sh and follows by interpolation arguments 
similar to those employed in the proof of [181 Lemma 5.2]: Define 0 = 1 — 2e (we as¬ 
sume e < 1/2). Then H-^ 2+e (n) = ( H 1 / 2 ~ £ (f2)y = ((L 2 (17), B^?({2)) g 2 )' = 
((-® 2 ^ 2 (^)) , j T 2 (l7)) 0 2 , so that the interpolation inequality yields \\v\\ H -i/ 2 + c ^ < 
The argument is completed by noting in view of suppw C 

Sh that 


IMU 2 (r 2 ) < IMU a (n) and \\v\\ {b i/ 2 W) , < CVh\\v\\ L 2^ n) , 

so that \\v\\ H -i/ 2 +^ n) < h (1 ~ e)/ -\\v\\ L 2 {s7) < h 1/2 ~ £ \\v\\ L 2 {n) . This proves (iBAl) . 

The bound (El follows from El and ||Au,— 77^ h Au,||^2( r ) < Ch Ikklir'^r) 
/i e ||w|| ff 3/2+e ( - r2 ) < h £ \\v\\ H -i/ 2 + C (ny An appeal to (IB.5D finishes the proof. □ 

We generalize □1 Lemma 5.5]: 

Lemma B.3 [generalizations of \18\ Lemma 5.5]] Assume that 17 satisfies Jf/Sj, 
Assumption (5.2)]. Then, for v G L 2 (Sh) C L 2 (l7) and w := T D (v) and the 
mortar approximation Wh of w, there holds 

IIV(it; - w^)||z,2 (r . L i ) < Ch 3/2 ( 1 + k,i| lnh|)||w||i,2(f2). 
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Proof Inspection of the proof of |fS Lemma 5.5] shows that we have to estimate 
the following two terms: 

Ti := / ||V(w - w h )\\ L *^) and t 2 ; = / IIV ( w — w ?i)I|l 2 (7 t) - 

J T=0 jT=ch 

Estimating Ti: Inspection of the proof of |18< Lemma 5.5] shows that 

Ti < Vh\\V(w - iy/0IU a (S( O ,c/O) < Vh\\V(w - w h )\\ L 2 (n .y 

We conclude together with (HOI) 

Ti < Ch 3/2 \\v\\ LHn) - 

We now turn to estimating T 2 . As in the proof of [18l Lemma 5.5], we consider 
the lowest order case k = 1 and the higher order cases k > 1 separately. 

Estimating T 2 for k = 1: Inspection of the proof of ns Lemma 5.5] gives 
(cf. [T8] eqn. (5.13)]) 

/ ||V(w — Wh)\\ L 2 ( lT ) dr (B.6) 

< |hi/i| 1/2 (||4 /2 V(w- llw)\\ L 2( C ,\ Scih) + ||^r 1/2 ( w - «h0llL2(C'VS cl fc))- 

The term w — I^w is estimated as in the proof of [18l, Lemma 5.5] with the 
aid of (weighted) Lt 2 -regularity asserted in [181 Lemma 5.4] and the estimate 
IMIs^ctt) < Vh\\v\\ L 2(Q) of [181 Lemma 5.2]. In total, we get 

- llw)\\ L 2 {C '\S aih ) < |ln/i| 1/2 /r 3/2 ||n|| L 2 (r2) . 

The second contribution of the right-hand side of (IbTgIi has to be treated with 
more care than in the proof of ns Lemma 5.5], where the convexity of 17 was 
conveniently exploited in order to control ||io — u>h\\L 2 ( 0 )'j more precisely, we do 
not have full 77 2 -regularity but only the limited shift theorem of [181 Assump- 
tion (5.2)]. In order to control the term ||<5 r ' (w — Wh)\\L 2 (C'\S clh ) appearing on 
the right-hand side of (IbUI) . we proceed by yet another duality argument. Let 
1 /) := T D (Sf 1 (w — Wh)), where 5r is the regularized distance function Sr ~ h + Sr- 
Note that Sr ~ Sr on C' \ S Cl h■ Denote by Ay, the Lagrange multiplier for ip, 
i.e., Ay,| 7l = — d n ip\n sm • From Lemma lB.41 we have for some s > 0 given by the 
stipulated shift theorem f |18l Assumption (5.2)]) 

IIV’I~ I ^h\ 1/2 \\S~ 1/2 (w - w h )\\ L 2 (n) , (B.7) 

+ ll^llff3/2 +e(r2) </i| £ \\5 r 1/2 (w-w h )\\ L * w . (B.8) 

The pair (ip, Ay,) solves the following saddle point problem: 

a(z,ip) + b(z, Ay,) = ( z,S~ 1 (w - w h )) Vz £ {z £ J^/f 1 (/7 i ) | = 0}, 

(B.9a) 
(B.9b) 


b(ip,q) = 0 Vq £ n^oo 2 ^)- 



30 


T. Horger et al. 


The Galerkin orthogonality satisfied by w — Wh reads for arbitrary lip G Vh 
a(w - Wh, lip) + b(Itp, X w ) = 0. 

Hence, we get by taking z = w — Wh in (IB.9al) 

\\Sp 1/2 (w - w h )\\L 2 (n) = a ( w ~ w h ,ip) + b(w - w h ,X^) (B.10) 

= a(w - Wh, ip ~ lip) ~ b(Iip, X w ) + b(w - w h , X^). 

We estimate |a(u> — Wh,ip — lip) \ < Ik — WhWw II'*/’ — IiPWh 1 - Since lip G Vh and 
[ip\ = 0: 

I b(Iip, Am)| = b(ip — lip, X w — n Mh Xw)\ < ||ip — Iip\\L 2 (r) || Atu — n^ h X w ||l 2 (.t) 

< Ik - m\%\n)U - ^Whi^WXwWh^d, 

where, in the last step, we employed the multiplicative trace inequality and the 
approximation property CEED- For the term b(w — Wh, A- 0 ), we employ again that 
[u>] = 0 and that Wh G Vh to get 

b(w-Wh,X^)= I [w - wh](X^ - IlM h X^) 

= I ([w - Wh] - IJM h [w - W h ])(X^ - IlM h X^) 

< h 1/2 \\[w - Wh]\\ H i/2^h e \\X^\\ H e^ r ). 

Noting ||[u; — W/jIH „i/ 2 , r , < Ww—WhWH 1 we obtain by inserting the above estimates 

in (IbTToIi 

\\S~ 1/2 (w - w h )\\ 2 L 2 (n) < 

ik - u> h \\m i \ip - iip\\m + ik - 

+ h 1/2 \\w - Wh\\mh e \\X^\\H^- 

The terms involving ip—lip are estimate with the aid of Lemma[B3]and the bound 
||V'|| s 3 / 2 ( fi ) given in (IB.71) : the terms ||w — wh\\hi are controlled in (IB. 31) : using 

ll A kk'„(r) ^ lk||H-i/a+«(n) and dEU) for ||A^|| H * a (r) we get 

II 5~ 1/2 (W - W h )\\ L 2 < 

h 3/2 \ ln/i| 1 / 2 ||u|| L 2 (r2) + h 1 ] ln/i| 1 / -||u|| L 2 (r 2 ) /i £ ||u|| ff -i/ 2 + e(fi) + h 3/2 h £ h~ z ||u|| L 2 (r2) . 
Now the result follows from (iBAl) . 

Estimating T 2 for k > 1: We proceed similarly to the case k = 1. The 
difference is that (cf. [181 (5.12)]) we need to control 

h- £ \\S- 1/2+e (w-W h )\\ L 2 {c , XSsih) . 

As in the case k = 1, we set up a dual problem with solution ip = T D (5p 1+2e (w — 
Wh)) and corresponding Lagrange multiplier ABy Lemma lB.41 (and trace esti¬ 
mates) we have the regularity assertions 

Il'Mk'^CO + Ilk'll HV 2 + t ((j) < ||5 r 1/2+£ (w - w h )\\ L 2 {n) . (B.ll) 
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Proceeding in the same manner as in the case k = 1, we arrive at 

||ir 1/2+e O - w h )f L ^ n} 

< \\w - w h \\ H i\\ip - IiJj\\ h i + H ^ I^\\Hih e \\Xw\\H‘ w (r) 

+ h 1/2 \\w - w h \\ H ih e \\\^\\ H e pw - 

The regularity assertions (IB. 1 ip as well as the approximation properties of Lemma|R5] 
yield ||<5 -1/2+£ 0 - w h )< h 3/2+e \\v\\ L 2 ^ n) and hence 

h~ z || S~ 1/2+£ (w - w h )\\ L 2 (c ,\ S; , ih) < h 3/2 \\v\\ L 2 (n) . 


□ 

Lemma B.4 (Generalization of [18L Lemma 5.2]) Let 5p be the regularized 
distance function. Then for the operator T D we have for w := T D (5pv) (and 


e > 0 sufficiently small): 

IMI B 3£ ( „) < C|ln/ l | 1 / 2 ||5- 1/2 n|| L2(f2) , (B.12) 

Hlfr*/*+* ( fl) < C e h- £ \\5f 1/2 v\\ L 2 (n) . (B.13) 

For w = T d (5f 1+2s v) we have 

IMIH3/ 2 +»(f2) < C\\5 r 1/2+£ v\\ L 2(n)- (B-14) 

Proof The proof is done with the same arguments as those of Lemma l2.8l □ 


We need an approximation result for the approximation from the constrained space 
Vh for functions that do not permit nodal interpolation. 

Lemma B.5 (approximation from constrained space) Let the constrained 
space Vh be defined in (18i (2.4b)] and assume hypotheses fltft (Al), (A2)J. Define 
the operator Ph as in \18[ ()-2)j but replace the interpolation operator 1^ by a (sub- 
domainwise) Scott-Zhang operator I z (cf. J2(fj ) that conforms to the boundaries 
of the subdomains fli, i = 1,... ,M. Then Ph is defined on H s (f2) PI Hq(T 2) for 
s > 1, it maps into Vh, and has the approximation properties 


^E ll u _ < ch s ^Miifs^), 


1 - Phv\\ L 2(n) < Ch s \\v\\ H s {n) , 1 < s < k + 1, 


E ii u - W) < ch a ~ r ii«ii j b 5i0o( « ) . 


Ik - Phv\\ L 2 {n) < Gh s |k|| B | i0o(fi) , 


1 < s < k + 1, 


s 0 N. 
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Proof We only show the first estimates as the estimates with B| j00 (f7)-regularity 
follow by interpolation arguments. It suffices to study the contribution Ekllh [I k ' SZ v\ 
to the operator P},. We recall that by the multiplicative trace inequality ||Hli 2 (ar 2 4 ) < 
ll u; l|L 2 (i7 i ) IIHI-Hpj?*) and the simultaneous approximation properties of I k ' SZ in 
L 2 and H 1 we have ||u — I k ' SZ < C'/i s_1/J ||u||H»(i 7 i )- Exploiting the 

Testability of the mortar projection 77^, we get the bound 

\\E k n h [I k ' sz v}\\ m(ni ) < h- 1 h 1/2 \\[I k ’ sz v]\\ LHani) 

<h 1 ■■||u - I k SZ v\\ L 2( a n i ) < h s 1 1|^|| H s {f2i )■ 


□ 
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C Details concerning Lemma 12.91 

The following lemma is an expanded version of Lemma [2.9l It is closely related to 
Lemma r2.3l with /3 = 1/2 there. The essential difference is that we replace the norm 
IMIff3/ 2 (f2) with the weaker norm ||2i|| B 3/2 ^ at the expense of a factor | In/r| 1 ^ 2 . 

Lemma C.l Let the bounded Lipschitz domain ft C R d , d G {2,3}. Assume that 
w G B 2 ^(ft) is a solution of 

—V ■ (AVw) = v. 

(i) There exist constants C (depending only on the lower bound oo for the eigen¬ 
values of A, the norm || A|| c , 0 ,i( 75 p and ft) and ci > 0 (depending only on 
ft) such that with the distance function Sr 

|| \ArV 2 w|| L 2 (r2 \ Sft) < C\J | ln/i||M| B 3/2 + C\\ v / 5r'f|U 2 (r2\s cl h)- (C- 1 ) 

In particular, ifv\ n \s aih = 0 then || \/frV 2 w\\ L 2(n\S h ) < C^/J\nh\\\w\\ B 3/^ ny 

(ii) Let c' > 0 be fixed. Assume t>|r 2 \s , = 0. Then there exist c, c 2 > 0 (depend¬ 
ing only on ft) such for every a > 0 

||<5?V 3 HU2 (fASa0 < C A , n [||5r 1 V 2 W || L 2 (fi \ Sc2h) (C.2) 

+ ll A llc 1 a(t2)||<5rVu;|| Z/ 2( fi \ Sc2h )j; 

the constant C A ,n depends on the coercivity constant Qo and || A|| c , 0 ,i( 7 j) as 
well as ft. 

(Hi) Let d > 0 be fixed. Assume , = 0. Assume that w G H' i ^ 2+e [ft) for 

some e G (0,1/2). Then there exists c > 0 (depending only on ft and c ) such 
that 

l|V 2 w|| L 2( fi \ Sa ) < C A , e h 1/2+e ||w|| if 3/2+ e(f 2), 
where C A , e depends on ao, 11A11c°,i( 7 ^)? and e. 

Proof 

Proof of & : We may restrict our attention to a local situation near a part of the 
boundary. The boundary r can locally be described by a graph </. That is, in a 
suitable coordinate system, we can define cylinders 

Cs = {(x,<j)( x) + t)\5 <t<D +5, x£ B}, 

C's = {(x, 4>{x) + t)\5 <t < D' + 5, x e B'}, 

where 5 < D < D' and B, B' are two concentric balls with B CC B'. Furthermore, 
for 5 = 0 we assume C'o C ft. In particular, {(x, <j>(x)) \ x G B'} C r. We also note 
that t ~ dist((x, 4>{x) + t), T). 

Let C's' be a third cylinder of the form C's = {(x, 4>(x) +t): 5 < t < D" + 5, x G 
B"} where B CC B" CC B' and D < D" < D'. Let y G C' 00 (M d ) be such that 
xlcf = 1 and x|r2\c(, = 0. To simplify the notation, we assume that the functions 
w, A, v are given in a coordinate system commensurate with the coordinate system 
describing the cylinders Cs, C's, viz., w evaluated at a point (x,<f>(x) + t) G C’ s is 
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given by w(x,4>(x) + t). A translation in the last variable defines the function w 
by ui(x, <t>(x ) + t ) := w(x , 4>(x) + t + 25). We note 

- V • (AVui) = v in C'_ 2 &\ (C.3) 

here (again in the coordinate system used to describe the cylinders) A{x,cj){x) + 
t ) = A(x , 4>{x) + t + 25) and v(x, <j>{x) + t) = v(x, 4>{x) + t + 25). 

1. step: We show (if 5 is sufficiently small) 

IMIh3/2 ( c") < Ca.o [v^MII^II^ (fl) + tf||®|U a (CL,)] , (C.4) 

where Ca.q depends on 17, the coercivity constant ao and II A|| c o,i( 72 )- 

Using the characterization of H :i ' 2 (Co) = {H 1 {Cq), i7 2 (Co))i/ 2 ,2 in terms of 
the A-functional, we write (cf. also |5l p.193, eqn. (7.4)]) 

IIxw||h 3/2 (c -) = j^(t~ 1/2 K(t,xw)) J 

= [ (t~ 1/2 K(t,xw)) Y + j (t~ 1/2 K(t,xw)) y-(C.5) 

•J t=0 J t = £ 

The second integral in (1C.51) can be estimated by 

[t~ 1/2 K{t,xw) S j Y ~ J t < ln e||xw||^3/2 (c ,). 

For the first integral in (1C.51) we employ interior regularity estimates for solutions 
of second order elliptic equation with vanishing right-hand side. Specifically, (ED 
and interior regularity (see, e.g., pi Thm. 8.8 and proof]) give (here, we assume 
that 5 is sufficiently small) 

llX^llff^C') < C A [<5 ^Iwllifl^j) + ||?||l, 2 (C'_ s )j ! 

where the constant Ca depends only on the coercivity constant ao and ||A|| c0 ,i( 7 j). 

Hence, estimating K (t, yw) = inf „ 6H 2 \\xti-v\\m( C ' 0 )+t\\v\\ H 2 ( C ' 0 ) < i||xfa||// 2 (C'), 
we obtain 

t~ 2 K 2 (t, x«5) dt < £\\xw\\ 2 H 2 {c - o) < C A \\w\\m(C'_ s ) + £ Ml^c_ s )] ■ 

We conclude 

!Ixw||h3/2 (C ') < c [£<W 2 ||w||?fi(c^) + e||v||| 2 (cli) + lne||x«j|| 2 3/2 (£,,)] (C.6) 

< Ca,si | “ + Inej ||w|l b^(c') + £ ll^ll?- 2 (cc 5 )} > 

where, in the last step we have employed that multiplication by a smooth function 
and translation are bounded operations on Sobolev (and therefore also Besov) 
spaces; the constant Ca,q depends on 17, ao, and ||A||c-o,!( 72 )- Selecting e = 5 2 

shows ||x6?||jt3/2 (C /) < C A ,n [\/l ln <5||MI B 3/2 (n) + ^NU 2 (cr 4 )] from which we 
get (1C.4I) in view of the support properties of x- 
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Step 2: Let z solve, for a fixed p > 0 and a parameter r < 1, the equation 

—V ■ ^AVz^j = / in a ball -B r (i+ P ) of radius r(l + p). 

We claim: For a constant C a, p > 0 that depends solely on the coercivity constant 
ao, IIA|| c o,i( 72 ); and P there holds 

I|V 2 z||l 2 ( Br ) < Ca,p [||/I|l 2 (b k1+p) ) + ||Vz||i 2 (B r(1+(>) ) + r ~\z\h*/ 2 (b T ( 1+p) )\ >( c - 7 ) 
||V 2 2 !||L 2 ( B r) < Ca, p [||/|U 2 (B r(1+p) ) +r- 1 ||V3|| i 2 ( B r(1+p)) ] ■ (C.8) 

We point out that the H 3 ' 2 -seminorm in is defined in terms of the Aronstein- 

Slobodeckij norm for scaling reasons. The bounds (1C.71) . (1C.81) follow from interior 
regularity in the following way. 

Scaling the ball 77 r (i+ P ) to a ball Bi+ P of radius 1 + p leads to an equation of 
the form 

— V • ^AVzj = r 2 / in a ball B\+ p of radius (1 + p), (C.9) 

where A and / are the coefficient and the right-hand side in the scaled variables. 
We note that 


l|V J A|| Loo(Si+p) ~ riV J A|| it » (BW , j e N 0 . 

Then standard interior regularity (see, e.g., [9j Thm. 8.8]) gives 



where the constant Ca, p depends only on the coercivitiy constant of A, the norm 
IIA||co,i(b 1+ ) j and p. In view of r < 1 and the fact that the operator A is ob¬ 
tained from A by an affine change of variables (a translation and an orthogonal 
transformation) we bound II A|lc?o,i(_gi + ) — ^ll^llc°’ 1 (72) an< ^ § e ^ 

ll ^ 7 Z IIl 2 (Bi) — Ca,p II/IIl 2 (B 1+p ) + II^IIh^Bj+p)] > (C-ll) 

for a constant Ca, p that depends only on p, the coercivity constant of A, and 
II A|| cr 0 , 1 ( 72 ) • Since the constant functions are in the kernel of the operator —V • 
(AV z) it is easy to conclude with a Poincare inequality that (1C.Ill) implies 

I|V 2 % 2(Si) < Ca, p [r 2 ||/|| L 2 ( B 1+p) + l|V% 2 (B 1+ J ■ (C.12) 

Scaling this equation back to S r ( 1+p ) yields the desired bound (1C.81) . For the proof 
of (1C.71) . we have to bring in the H 3 ' 2 -seminorm. Let n G Vi be arbitrary. Then 
the function z — tv satisfies in view of the fact that X7n is constant 

-V- (AV(z —tt)) = r 2 / + V • (AW) =r 2 / + (V-A)-V7r=:/. 
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(We employed the convention that the divergence operator V- in the expression 
V • A acts on columns of A). Applying (1C.1111 to this equation (and replacing r 2 / 
with / and z with z — tt) yields 

ll^ 7 Z W i 2 (Bi) — W Z ~ 7F II [|I/IIl 2 (B 1+p ) + W Z ~ 7r II.W 1 (B 1+p )] 

< Ca,p [r 2 \\f\\LHB^ 1+p) ) + r\\Vn\\ L2{Si+p) + \\z - tv\\h 3/ 2 (Si+p) \ 

< Ca,p [r 2 \\f\\L*(B H 1 +p) )+r\\Vz\\ LHSi+p) + \\z-tt \\ h3/a(gi+p) ] . 

Infimizing over all n G V\ yields 

ll^ 7 z \\ B 2 (Bi) — Ca,p [p II/IIl 2 (B 1+p ) + r ll^ 7 ^lll, 2 (B 1+p ) + Mb 3 / 2 (B 1+p )] ■ (C.13) 

Scaling back to B r ( 1+p ) yields (1C.71) . if we note the scaling properties of the 
Aronstein—Slobodeckij seminorm. 

3. step: Applying the result of step 2 to the function w yields 

l|V 2 w|U 2 (B r ) < C A ,p [l|Vw|| L 2 (B r(1+p) ) +r 1/2 \w\ H 3/2( Ba+p)r) + ||?||l 2 (b p(1+p) )] 

(C.14) 

for all balls B r such that B^ 1+p y C C'_ ag . Using, for example, the Besicovitch cov¬ 
ering theorem [3, we can cover Co by overlapping balls B ri (xi) with centers Xi and 
radii r,; ~ 5r(xi) such that the stretched balls B r .^i +p ^(xi) have a finite overlap 
property (see [181 Lemma A.l] for details). A covering argument and afterwards 
(ICAll then show 

II \AtV 2 ui|| L 2( Co ) < C A ,n [|MIb 3/2 (C //) + || \ArVw|| l 2 (c') + II \f&rv ||l 2 (c')] 

< C A ,n [a/| ln£||M| B 3/2 (n) + s \\v\\l 2 (C’_ s ) + II v/rv|| i 2 (C / i) j , 

where the constant Ca,q depends only on the coercivity constant ao, ||A|| c0 ,i( 7 jp 
and fl. 

Since w and v are obtained by a translation, we arrive at 

II\/^V 2 W |U 2 (C2S) < C A ,n [v / |ln5|||w|| B 3/2 (r2 > + || \Ar«IU 2 (fi\Sc 2S )] 

for a suitable C 2 that depends solely on the Lipschitz character of r. Taking <5 ~ h 
produces the desired result. 

Proof of {§[): The estimate merely expresses interior regularity for solutions 
of elliptic equations with vanishing right-hand side. It follows from (1C.81) and a 
covering argument. More precisely, as in Step 2, we start from 

—V • (AVz) = f in a ball B r (i+ p ) of radius r( 1 + p). 

Differentiating this equation once gives for a G with |a| = 1 

—V ■ (A VD a z) = D a f + V • (D a A) Vz in a ball B^+p) of radius r( 1 + p), 
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where again the divergence operator in the expression V • (D a A) acts on the 
columns of D a A. We get from (1C.81) by considering all a G N d with |a| = 1 

l|V 3 z||z,2( Br .) < 

Ca, p [||V/|U 2(Br(1+p)) +r- 1 ||V 2 ^|| i 2 (Br(i+p)) + ||V 2 A|| B oo (Br(1+p)) ||Vz|| B 2 (Br(1+p)) j . 

A covering argument then produces the claim since / is assumed to satisfy / = 0 
on 12\ S c ’h in the statement Jn]). We also note that || V 2 A|| B oo(f 2 ) can be bounded 
by || VA|| w i, «,(„). 

Proof of ifrnl) : Our starting point is the above Step 2: We claim that for the 
function z satisfying (IC.9I1 we have 

l|V 2 2||z,2 (Br ,) < C A [||/||l 2( Bt . (1+p)) + ||Vz|| L 2 ( Br(1+rt) +C 1/2 + E |2| B 3/2 + e(Br(1 + p)) j , 

(C.15) 

where the H 3 ' 2+e -seminorm is again an Aronstein-Slobodeckij norm, and the con¬ 
stant Ca depends only on the coercivity constant ao, ||A|| c 0 ,i( 77 p and e G ( 0 , 1 / 2 ). 

To see this, we check the derivation of (1C. 131) . We see that one can (marginally) 
modify the arguments to obtain instead of (1C. 1311 the estimate 


|V 2 ?| 


L 2 (B 1+p ) 


< 


2 (B 1+p ) + r ll Vz ll L 2 (B 1+p ) + \ Z \hW+c(B 1+p )’ 


where the implied constant depends on p, ao, 11A||co,i( 77 ), and e. Scaling back to 
B r ( i+p) yields 


r 2 ||V 2 2 || L 2 (Br(1+p)) <r 2 ||/|| Zl 2 (Br(i+p)) +rr||Vi : || B 2 (Br , (1+p)) +r 3/2+£ |^| ff 3/2 +e(Br(i+p)) , 

which leads to (1C. 1511 . We now use / = 0 on (2 \ S c 'h and use these estimates with 
w in place of 2 . A covering argument then gives 


II V 2 wl 


M lk 2 (f2\S ?fc ) ^ ||Vu;|| L 2 ( fi \ Sc , h ) + h 1/2+e ||u;|| B 3/2+ e ( f2) . 

This concludes the proof of (Hull . 
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D Details of [24, Sec. 5.3] 

By and large, we follow the arguments of [241 Sec. 5.3] and adapt them as needed. 

In what follows, T is a quasi-uniform mesh with mesh size h. We consider the 
bilinear form 



where the matrix A is sufficiently smooth and pointwise symmetric positive defi¬ 
nite. 

D.l The interior case 

In this subsection, we assume that Vh = S k,1 (T). 

For balls Bd of radius d, it will be convenient to introduce the notation 


IMIkfid M Hi(B d ) + ^IM \L 2 (B d )- 


We start with making the notion of “superapproximation” more precise: 

Lemma D.l (superapproximation on balls) Let Bd C 17 be a ball of radius 
d. Let uj G C°° with suppo; C f?d /2 an d 


liv'wiu- < Cd~\ J = 0,...,k. 

Then for every u G 14, the interpolant I(uj 2 u) G S k ’ 1 [T) fl Hq(Ba) satisfies 
\u 2 u- I(u 2 u)\ m{Bd) < C-\\\u\\\ 1}Bd , 


(D.l) 


(D.2) 



(D.3) 


We assume implicitly that d > h is sufficiently large. 

Proof Since suppw C B ^/2 and d is large compared to h, we have supp/(w 2 u) C 
Bd- For each element K we have the estimates 


\\uj 2 u - I(uj 2 u)\\ L 2( K ) + h\\X7(uj 2 u - I(uj 2 u))\\ L 2( K ) < C h k+ 1 \\\7 k+1 (w 2 u)\\ L 2 (k)■ 


Inductively, we see that ||V : '(a; 2 )||z / oo < Cd J for j = 0,..., k. Using the fact that 
u is piecewise polynomial of degree k we conclude 


k 


|| (X7 k+1 (ui 2 u)\\l2(k) < C^2d~ (k+1 ~ j) \\V j u\\ LHK) . 


3=0 


An inverse estimate produces in view of h/d< 1 



□ 
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The next result shows an inverse estimate for “discrete harmonic” functions: 

Lemma D.2 Let uu G Vh satisfy 

a(uh,v) = 0 Vv G Vh with suppt C (D.4) 

for a ball Bd C fl of radius d (implicitly assumed sufficiently large compared to 
h). Then 


\\Vu h \\ LHBd/2) < Cd 1 \\u h \\ L i( Bd y 

Proof Select a cut-off function uj with supp w C Bd and da as well as x = 1 on 
Bd/ 2- 

We write for arbitrary x € Vh with supp x C Bd 


f 

J n 


cj AVuh • Vu/i = / A Vuh 


— j 2uf l Lj'V(jj • A S7uh 
Q 


/ AViijj • V(cu 2 Uh) — / 

J £2 J,ft 

= / AViih • V(oj 2 m/j — x) — / 2uhU}\7u) ■ ATJuh- 

■In Jn 


We conclude from Lemma lD.il with x = I{ui 2 Uh) there and Young’s inequality 


/ 


uj 2 AVu h ■ Vu h < ^||Viih||i, 2 (Bil) ||u fc ||| 1)BiJ + 


We conclude 


l|v^|| i2(Bd/2 > < y ^|||uh|||i ,Bd + -f\\ u h\\L 2 (B d )- 

Iterating the argument yields 

||VU/ 1 || Z/ 2( Bd/4 ) < -||w^|||l,B d + -j\\uh\\L^(B d )- 
Finally, a standard inverse estimate produces 

||Vu/i||i,2( Bd/4 ) ^ ~\\ u h\\L 2 (B d ), 

which is the desired final bound. □ 

We now show the main result: 

Theorem D.3 Let u G H 1 (f2) and uh € Vh be such that 

a(u — Uh,v) = 0 Vv G Vh with suppn C Bd 

for a ball Bd C 17 of radius d (implicitly assumed sufficiently large compared to 
h). Then 

\\V{u-u h )\\ L 2 ( B )<C inf lllu-xlli.Bd +Cd~ l \\u- Uh\\ L *(B d ) (D.5) 

XGVh 

<C inf ||V(u-x)||z, 2 (B d ) +C'cT 1 ||u-u^|| L 2 (Sd) . (D.6) 
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Proof The second bound (ID. 61) follows from (ID. 51) by an application of the (second) 
Poincare inequality. 

Let La G C°° with suppy C B d be such that w = 1 on B d / 2 and assume (ID. II) . 
Define 

u := cau 

and let Uh G Vh fl Ho(B d ) be its Galerkin approximation on Vh fl -Bo(7?d), i.e., 
a(u-Uh,v) = 0 Vv G Vh D Ho(B d ). 

Then it is classical (and easy to see) that 

/ V5/i • (AVi/,) < / Vu- (AV«). 

JB d JB d 

In particular, we get 

\uh\m(B d ) < C\u\ H i( Bd ) < |||u|||i,B d . (D.7) 

Next, we write 

u- u h = (u-u h ) + (u h - u h ) in B d /2 (D.8) 

and estimate each of the two terms separately. For the first one, we employ (ID. 71) 
to get 

I|V(m- u h )\\ L 2(B d/2 ) ^ ||Vu|| L 2 (Bd) < |||u|||i, Sd . (D.9) 

For the second term in (ID. 81) . we observe that Uh — Uh is discrete harmonic in B d / 2 
(in fact, almost in B d ) since 

a(uh ~ Uh,v ) = a(u - u,v) = 0 Vv G Vh D Ho(B d / 2 ). 

Therefore, Lemma ID. 21 is applicable and yields in view of ui = 1 on B d / 2 : 

l|V(5 h - Uh)\\L 2 (B d/4 ) < Cd 1 \\u h - U h || L 2 (S d/2 ) 

< Cd~ 1 \\u h - u\\ L 2 ( Bd/2 ) + Cd^Wu - u h \\ L 2 ( Bd/2 ) 

= Cd~ 1 \\uh - u\\L*(B d/2 ) + Cd^Wu - u h \\ L 2 (Bd/2) , 

where in the last step we exploited u|s d/2 = u\ Bd/2 due to cu| B<J/2 = 1. Since Uh, 
u G Hn(B d ), a Poincare inequality together with (iDTl) produces 

II V(Wfc - U h )\\ L 2(B d/i ) < ||V(Uft, - u)|| L 2( Sd ) + d 1 ||l* — Uh || L 2 (B d / 2 ) 

^ \\\u\\\i,B d +d~ 1 \\u-u h \\ L ^ Bd/2) (D.10) 

Combining (E3 and (ID.1011 yields again with ca = 1 on B d / 2 

l|V(« - Uh) |U=(S d/4 ) < l|V(u — Uh)\\ L 2 (B d/i ) + ||V(u h - U h )\\ L 2(B d/i ) (D.ll) 

< IMI|l,-Bi + d 1 H U “ U h \\ L 2( Bd ). 

The final step consists in noting for arbitrary x G Vh that u—uh = (u— x) + (x — u h) 
so that an application of (ID. Ill) applied to u — x yields 

IIv(it - U h )\\ L 2(B d/i ) ^ |||w - xllll, B d +d 1 ||u U h \\ L 2(B d )- 

□ 
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D.2 The boundary case 

We now check to what extent the above results extend up to the boundary. We 
Let r be (a part of) the boundary dfi. We denote by 

Dd = BdC\ f2 

semiballs of radius d, where the implict assumption is always that the center of 
Dd lies on r. Another implicit assumption in the following is that for all semiball 
appearing in the following, we assume 

dD d ndncr 

We define |||'rr|||: = Mff 1 (D d ) + l\\ u \ U 2 (D d )- We employ a space Vh. which sat¬ 
isfies 

Vh C S k,1 (T), u| r = 0 VveV h . (D.12) 

We employ the following observation: For a semiball D d C B d (same center) as 
described above, oj G C(f(Bd), and u G Vh, we have urn G Ho(Dd). 

Compared to Lemma fD.41 the cut-off function lj may be = 1 near parts of r: 

Lemma D.4 (superapproximation on semiballs) Assume W.12\) . Let D d C 

fi be a semiball of radius d. Let ui G C°° (M ra ) with supp ui C Bd /2 and W.l\) . Then 
for every u G Vh the interpolant I(ui 2 u ) G S k ’ 1 ('T) fl ILo(Dd ) satisfies 

\u 2 u - I{lo 2 u) |m(D d ) < C^IMIIi.A*, (D.13) 

i||w 2 M- I(u 2 u)\\ L 2 {Dd) < C ||w|||i,D d . (D.14) 

We assume implicitly that d > h is sufficiently large. 

Proof Follows by the same arguments as in Lemma ID. II □ 

The next result shows an inverse estimate for “discrete harmonic” functions: 
Lemma D.5 Assume \D.ltA) . Let Uh G Vh satisfy 

a(uh,v) = 0 Vv G Vh with supp PC Dd (D.15) 

for a semiball Dd C J? of radius d (implicitly assumed sufficiently large compared 
to h). Then 

l|Vu/ l || Z; 2( £ , (i/2 ) < Cd 1 ||u^|| Z/ 2( £ , d ). 

Proof Again, this follows by tracing the arguments in the proof of Lemma ID. 5 1 □ 

Theorem D.6 Assume \D.12\) . Let u G Pd 1 ^) with u\r = 0 and Uh G Vh be 
such that 

a{u — Uh,v) = 0 Vu G Vh with suppw C D d 
for a semiball Dd C f2 of radius d (implicitly assumed sufficiently large compared 
to h). Then 

l|V(u — Uh)\\ L 2 (D d/4 ) < C infk \\u-x\\i,D d +Cd~ 1 \\u-Uh\\L^D d ) 

< C ll V ( M - X)h^D d ) + Cd-'Wu - U h \\ L ^D d )- 

Proof The second inequality follows again from a Poincare inequality (the first 
one, this time). The first inequality follows again from tracing the arguments of 
the proof of Theorem ID. 31 □ 
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E Detailed numerics 

E.l Slit domain 

The geometry is a slit domain 

C=C S :=(-l,l) 2 \[0,l)x{0}. 

The exact solution is given by 

\x — Xo\ a sm(an(f>) 

for different choices of the parameters a and Xo (and a). 

The inhomogeneous Dirichlet boundary conditions are realized by nodal inter¬ 
polation. The equation considered is 

-Au = f. 

Starting from a coarse mesh, we perform a sequence of uniform (red) refinements. 
We consider a lowest order discretization, i.e., k = 1. 

Strictly speaking, the slit domain is not covered by our theory. Also not cov¬ 
ered by our theory are the variational crimes associated with approximating the 
inhomogeneous Dirichlet data. Nevertheless, we expect the convergence behavior 
detailed in Corollary 14.31 to be an good description of the actual convergence be¬ 
havior. We assume that the global regularity of the solution u is described by 
s = 1 + a (actually, it is 1 + a — £ for all e > 0). Corollary 14.31 then lets us expect 
for the two cases xo = (0,0) and xo ^ (0, 0) the following convergence rates: 

Xo = (0,0) => r = min{2,1 + a, — 1 + 1/2 + (1 + a)} = min{2,1/2 + a} 

Xq (0,0) => r = min{2,1 + a, —1 + 1/2 + oo} = min{2,1 + a} 

In the following Tables lETllE. 71 we vary the parameter a. In each table sepa¬ 
rately we vary the location. The locations under investigation are ( xo , yo) = (0, 0), 
(xo,yo) = (0.5,0) and ( xo,yo ) = (0,1). We observe that the theoretical conver¬ 
gence rates are mostly achieved in our numerical simulations. 



(%o,yo) = ( 0 , 0 ) 

a = 7r/2 

Oro, Vo) = (0.5,0) 
a = 7r 

{x 0 ,yo) = (0,1) 

a = 7r 

DOFs 

L 2 -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

2.9124e-02 

- 

3.8405e-02 

- 

3.0468e-02 

- 

348 

1.5745e-02 

0.89 

1.0451e-02 

1.88 

1.2883e-02 

1.24 

1.315 

8.1422e-03 

0.95 

4.8926e-03 

1.10 

5.2831e-03 

1.29 

5.109 

4.1322e-03 

0.98 

2.1508e-03 

1.19 

2.0814e-03 

1.34 

20.137 

2.0799e-03 

0.99 

8.1046e-04 

1.41 

7.9896e-04 

1.38 

79.953 

1.0430e-03 

1.00 

3.0969e-04 

1.39 

3.0187e-04 

1.40 

318.625 

5.2221e-04 

1.00 

1.1780e-04 

1.39 

1.1288e-04 

1.42 

1.272.129 

2.6125e-04 

1.00 

4.1750e-05 

1.50 

4.1903e-05 

1.43 

5.083.777 

1.3066e-04 

1.00 

1.4985e-05 

1.48 

1.5472e-05 

1.44 


Table E.l Slit domain, k = 1: Influence of the position of singularity for a = 0.5. 
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(xo,yo) = (0,0) 

a = 7r/2 

(xo,yo) = (0.5,0) 
a — 7r 

{xo,yo) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

1.1642e-02 

- 

1.6577e-02 

- 

1.6031e-02 

- 

348 

5.6371e-03 

1.05 

6.5153e-03 

1.35 

6.0796e-03 

1.40 

1.315 

2.6059e-03 

1.11 

2.3569e-03 

1.47 

2.2272e-03 

1.45 

5.109 

1.1835e-03 

1.14 

8.4203e-04 

1.48 

7.8285e-04 

1.51 

20.137 

5.3296e-04 

1.15 

2.7392e-04 

1.62 

2.6776e-04 

1.55 

79.953 

2.3888e-04 

1.16 

9.1902e-05 

1.58 

9.0042e-05 

1.57 

318.625 

1.0679e-04 

1.16 

3.0823e-05 

1.58 

2.9942e-05 

1.59 

1.272.129 

4.7666e-05 

1.16 

9.6773e-06 

1.67 

9.8788e-06 

1.60 

5.083.777 

2.1257e-05 

1.16 

3.1032e-06 

1.64 

3.2407e-06 

1.61 


Table E.2 Slit domain, k = 1: Influence of the position of singularity for a = 2/3. 



(xo,yo) = (0,0) 

a = 7t/2 

(x 0 ,y 0 ) = (0.5,0) 
a = 7r 

C xo,yo ) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

6.1391e-03 

- 

1.1088e-02 

- 

1.0692e-02 

- 

348 

2.8187e-03 

1.12 

4.1329e-03 

1.42 

3.8553e-03 

1.47 

1.315 

1.2351e-03 

1.19 

1.4164e-03 

1.54 

1.3388e-03 

1.53 

5.109 

5.3338e-04 

1.21 

4.7830e-04 

1.57 

4.4562e-04 

1.59 

20.137 

2.2846e-04 

1.22 

1.4725e-04 

1.70 

1.4420e-04 

1.63 

79.953 

9.7267e-05 

1.23 

4.6683e-05 

1.66 

4.5843e-05 

1.65 

318.625 

4.1233e-05 

1.24 

1.4761e-05 

1.66 

1.4401e-05 

1.67 

1.272.129 

1.7428e-05 

1.24 

4.3773e-06 

1.75 

4.4861e-06 

1.68 

5.083.777 

7.3524e-06 

1.25 

1.3285e-06 

1.72 

1.3889e-06 

1.69 


Table E.3 Slit domain, k = 1: Influence of the position of singularity for a = 0.75. 



{xo,yo ) = (0,0) 
a = 7r/2 

{xo,yo) = (0.5,0) 
a = 7r 

(xo,yo) = (0, l) 

a = 7r 

DOFs 

L z -e rror 

rate 

L z -error 

rate 

L z -error 

rate 

97 

4.1949e-03 

- 

3.0111e-04 

- 

2.9784e-04 

- 

348 

1.4605e-03 

1.52 

9.9257e-05 

1.60 

9.3618e-05 

1.67 

1.315 

4.8756e-04 

1.58 

2.9679e-05 

1.74 

2.8033e-05 

1.74 

5.109 

1.5909e-04 

1.62 

8.6201e-06 

1.78 

8.0205e-06 

1.81 

20.137 

5.1667e-05 

1.62 

2.2994e-06 

1.91 

2.2266e-06 

1.85 

79.953 

1.6874e-05 

1.61 

6.2533e-07 

1.88 

6.0606e-07 

1.88 

318.625 

5.5687e-06 

1.60 

1.6832e-07 

1.89 

1.6270e-07 

1.90 

1.272.129 

1.8596e-06 

1.58 

4.3035e-08 

1.97 

4.3240e-08 

1.91 

5.083.777 

6.2798e-07 

1.57 

1.1235e-08 

1.94 

1.1403e-08 

1.92 


Table E.4 Slit domain, k = 1: Influence of the position of singularity for a = 1.01. 
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(xo,yo) = (0,0) 

a = 7r/2 

(xo,yo) = (0.5,0) 
a — 7r 

{xo,yo) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

5.7534e-03 

- 

2.8888e-03 

- 

2.8799e-03 

- 

348 

1.9412e-03 

1.57 

9.1606e-04 

1.66 

8.6618e-04 

1.73 

1.315 

6.2583e-04 

1.63 

2.6267e-04 

1.80 

2.4731e-04 

1.81 

5.109 

1.9689e-04 

1.67 

7.2833e-05 

1.85 

6.7463e-05 

1.87 

20.137 

6.1446e-05 

1.68 

1.8673e-05 

1.96 

1.7871e-05 

1.92 

79.953 

1.9191e-05 

1.68 

4.8621e-06 

1.94 

4.6455e-06 

1.94 

318.625 

6.0229e-06 

1.67 

1.2518e-06 

1.96 

1.1921e-06 

1.96 

1.272.129 

1.9023e-06 

1.66 

3.0913e-07 

2.02 

3.0311e-07 

1.98 

5.083.777 

6.0474e-07 

1.65 

7.7720e-08 

1.99 

7.6552e-08 

1.99 


Table E.5 Slit domain, k = 1: Influence of the position of singularity for a = 10/9. 



(xo,yo) = ( 0 , 0 ) 

a = 7r/2 

(x 0 ,y 0 ) = (0.5,0) 
a = 7r 

C xo,yo ) = (o, i) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

97 

7.3549e-03 

- 

6.3401e-03 

- 

6.3849e-03 

- 

348 

2.2414e-03 

1.71 

1.8792e-03 

1.75 

1.7790e-03 

1.84 

1.315 

6.4849e-04 

1.79 

5.0365e-04 

1.90 

4.6905e-04 

1.92 

5.109 

1.8251e-04 

1.83 

1.3007e-04 

1.95 

1.1878e-04 

1.98 

20.137 

5.0718e-05 

1.85 

3.1798e-05 

2.03 

2.9443e-05 

2.01 

79.953 

1.4021e-05 

1.85 

7.8665e-06 

2.02 

7.2227e-06 

2.03 

318.625 

3.8699e-06 

1.86 

1.9356e-06 

2.02 

1.7642e-06 

2.03 

1.272.129 

1.0682e-06 

1.86 

4.6924e-07 

2.04 

4.3055e-07 

2.03 

5.083.777 

2.9514e-07 

1.86 

1.1524e-07 

2.03 

1.0519e-07 

2.03 


Table E.6 Slit domain, k = 1: Influence of the position of singularity for a = 4/3. 



(xq,Vo) = (0,0) 

a = 7r/2 

(xo,yo) = (0.5,0) 
a = 7r 

C xo,yo ) = (0,1) 
a = 7r 

DOFs 

L z -e rror 

rate 

L z -error 

rate 

L z -error 

rate 

97 

7.5901e-03 

- 

7.6006e-03 

- 

7.6553e-03 

- 

348 

2.1664e-03 

1.81 

2.1751e-03 

1.81 

2.0530e-03 

1.90 

1.315 

5.8638e-04 

1.89 

5.6614e-04 

1.94 

5.2238e-04 

1.97 

5.109 

1.5450e-04 

1.92 

1.4246e-04 

1.99 

1.2877e-04 

2.02 

20.137 

4.0197e-05 

1.94 

3.4615e-05 

2.04 

3.1388e-05 

2.04 

79.953 

1.0396e-05 

1.95 

8.5086e-06 

2.02 

7.6403e-06 

2.04 

318.625 

2.6803e-06 

1.96 

2.0921e-06 

2.02 

1.8651e-06 

2.03 

1.272.129 

6.8978e-07 

1.96 

5.1307e-07 

2.03 

4.5726e-07 

2.03 

5.083.777 

1.7730e-07 

1.96 

1.2691e-07 

2.02 

1.1258e-07 

2.02 


Table E.T Slit domain, k = 1: Influence of the position of singularity for a = 1.5. 
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E.2 L-shaped domain 


The geometry is an L-shaped domain: 

n mJi L := (— 1, l) 2 \ [0,1) x (-1,0]. 

The exact solution is given by 

\x — xo\ a sin(a 7 r/>) 

for different choices of the parameters a , xo, and a. 


E.2.1 Lowest order discretization k = 1 

Structurally, the situation is similar to the situation in Section IE. II From Corol¬ 
lary 231 we expect the following convergence rates: 

xo = (0,0) => r = min{2,1 + a, — 1 + 2/3 + (1 + a)} = min{2, 2/3 + a} 
xo / (0,0) => r = min{2,1 +a, —1 + 2/3 + oc} = min{2,1 + a) 

Also in this case the numerical rates depicted in Tables IE.8IIK. 1 3l are very close to 
the rates expected by our theory. 



(%o,yo) = ( 0 , 0 ) 

a = 7r/2 

(so, 2/0) = (0.5,0) 
a — 7r 

(xo,yo) = (0,1) 

a = tv 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

81 

1.1719e-02 

- 

1.0132e-02 

- 

9.5383e-03 

- 

289 

5.4059e-03 

1.12 

4.0619e-03 

1.32 

3.4014e-03 

1.49 

1.089 

2.3165e-03 

1.22 

1.4839e-03 

1.45 

1.1426e-03 

1.57 

4.225 

9.5790e-04 

1.27 

5.1844e-04 

1.52 

3.7569e-04 

1.60 

16.641 

3.8922e-04 

1.30 

1.7681e-04 

1.55 

1.2222e-04 

1.62 

66.049 

1.5663e-04 

1.31 

5.9408e-05 

1.57 

3.9513e-05 

1.63 

263.169 

6.2682e-05 

1.32 

1.9762e-05 

1.59 

1.2720e-05 

1.64 

1.050.625 

2.5002e-05 

1.33 

6.5263e-06 

1.60 

4.0824e-06 

1.64 

4.198.401 

9.9525e-06 

1.33 

2.1437e-06 

1.61 

1.3072e-06 

1.64 


Table E.8 L-domain, k = 1: Influence of the position of singularity for a = 2/3. 
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(xo,yo) = (0,0) 

a = 7r/2 

(xo,yo) = (0.5,0) 
a — 7r 

{x 0 ,yo) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

81 

6.1585e-03 

- 

6.8141e-03 

- 

6.2506e-03 

- 

289 

2.6986e-03 

1.19 

2.5648e-03 

1.41 

2.1211e-03 

1.56 

1.089 

1.1123e-03 

1.28 

8.8428e-04 

1.54 

6.7413e-04 

1.65 

4.225 

4.4037e-04 

1.34 

2.9202e-04 

1.60 

2.0903e-04 

1.69 

16.641 

1.7107e-04 

1.36 

9.4164e-05 

1.63 

6.4027e-05 

1.71 

66.049 

6.5689e-05 

1.38 

2.9909e-05 

1.65 

1.9471e-05 

1.72 

263.169 

2.5030e-05 

1.39 

9.4012e-06 

1.67 

5.8930e-06 

1.72 

1.050.625 

9.4877e-06 

1.40 

2.9328e-06 

1.68 

1.7774e-06 

1.73 

4.198.401 

3.5834e-06 

1.40 

9.0968e-07 

1.69 

5.3475e-07 

1.73 


Table E.9 L-domain, k = 1: Influence of the position of singularity for a = 0.75. 



(xo,yo) = (0,0) 

a = 7r/2 

(x 0 ,y 0 ) = (0.5,0) 
a = 7r 

(xo,yo ) = (0, l) 

a = 7r 

DOFs 

L z -e rror 

rate 

L z -error 

rate 

L z -error 

rate 

81 

4.6216e-03 

- 

1.8387e-04 

- 

1.6841e-04 

- 

289 

1.6860e-03 

1.45 

6.0370e-05 

1.61 

5.0364e-05 

1.74 

1.089 

5.4867e-04 

1.62 

1.8034e-05 

1.74 

1.3883e-05 

1.86 

4.225 

1.7284e-04 

1.67 

5.1378e-06 

1.81 

3.6942e-06 

1.91 

16.641 

5.2963e-05 

1.71 

1.4253e-06 

1.85 

9.6399e-07 

1.94 

66.049 

1.5970e-05 

1.73 

3.8870e-07 

1.87 

2.4842e-07 

1.96 

263.169 

4.7758e-06 

1.74 

1.0474e-07 

1.89 

6.3437e-08 

1.97 

1.050.625 

1.4238e-06 

1.75 

2.7971e-08 

1.90 

1.6086e-08 

1.98 

4.198.401 

4.2471e-07 

1.75 

7.4181e-09 

1.91 

4.0561e-09 

1.99 


Table E.10 L-domain, k = 1: Influence of the position of singularity for a = 1.01. 



(xo,yo) = (0,0) 

a = 7r/2 

(xo,yo) = (0.5,0) 
a = 7T 

(xo,yo) = (0,1) 

a = 7r 

DOFs 

L z -e rror 

rate 

L z -error 

rate 

L z -error 

rate 

81 

6.5660e-03 

- 

1.7641e-03 

- 

1.6229e-03 

- 

289 

2.3309e-03 

1.49 

5.5465e-04 

1.67 

4.6837e-04 

1.79 

1.089 

7.3413e-04 

1.67 

1.5847e-04 

1.81 

1.2424e-04 

1.91 

4.225 

2.2257e-04 

1.72 

4.3172e-05 

1.88 

3.1761e-05 

1.97 

16.641 

6.5650e-05 

1.76 

1.1458e-05 

1.91 

7.9588e-06 

2.00 

66.049 

1.9056e-05 

1.78 

2.9916e-06 

1.94 

1.9695e-06 

2.01 

263.169 

5.4810e-06 

1.80 

7.7249e-07 

1.95 

4.8311e-07 

2.03 

1.050.625 

1.5690e-06 

1.80 

1.9789e-07 

1.96 

1.1771e-07 

2.04 

4.198.401 

4.4822e-07 

1.81 

5.0393e-08 

1.97 

2.8532e-08 

2.04 


Table E.ll L-domain, k = 1: Influence of the position of singularity for a = 10/9. 
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(xo,yo) = (0,0) 

a = 7r/2 

(xo,yo) = (0.5,0) 
a — 7r 

{x 0 ,yo) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

81 

8.6776e-03 

- 

3.8962e-03 

- 

3.6446e-03 

- 

289 

2.8523e-03 

1.61 

1.1374e-03 

1.78 

1.0008e-03 

1.86 

1.089 

8.2870e-04 

1.78 

3.0272e-04 

1.91 

2.5331e-04 

1.98 

4.225 

2.3073e-04 

1.84 

7.7239e-05 

1.97 

6.2153e-05 

2.03 

16.641 

6.2539e-05 

1.88 

1.9331e-05 

2.00 

1.5073e-05 

2.04 

66.049 

1.6688e-05 

1.91 

4.7956e-06 

2.01 

3.6440e-06 

2.05 

263.169 

4.4099e-06 

1.92 

1.1852e-06 

2.02 

8.8167e-07 

2.05 

1.050.625 

1.1580e-06 

1.93 

2.9260e-07 

2.02 

2.1389e-07 

2.04 

4.198.401 

3.0279e-07 

1.94 

7.2263e-08 

2.02 

5.2069e-08 

2.04 


Table E.12 L-domain, k = 1: Influence of the position of singularity for a = 4/3. 



( xo,Vo ) = (0,0) 

a = 7r/2 

(x 0 ,y 0 ) = (0.5,0) 
a = 7r 

C xo,yo ) = (0,1) 

a = 7r 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

81 

8.9932e-03 

- 

4.7178e-03 

- 

4.4942e-03 

- 

289 

2.8151e-03 

1.68 

1.3287e-03 

1.83 

1.2166e-03 

1.89 

1.089 

7.8034e-04 

1.85 

3.4367e-04 

1.95 

3.0580e-04 

1.99 

4.225 

2.0751e-04 

1.91 

8.5903e-05 

2.00 

7.5035e-05 

2.03 

16.641 

5.3910e-05 

1.94 

2.1227e-05 

2.02 

1.8321e-05 

2.03 

66.049 

1.3835e-05 

1.96 

5.2331e-06 

2.02 

4.4827e-06 

2.03 

263.169 

3.5256e-06 

1.97 

1.2917e-06 

2.02 

1.1011e-06 

2.03 

1.050.625 

8.9467e-07 

1.98 

3.1955e-07 

2.02 

2.7158e-07 

2.02 

4.198.401 

2.2641e-07 

1.98 

7.9238e-08 

2.01 

6.7212e-08 

2.01 


Table E.13 L-domain, k = 1: Influence of the position of singularity for a = 1.5. 
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E.2.2 Second, order discretization k = 2 

All calculations are performed for xo = (0,0) and a = 2/3. Only the singularity 
parameter a is varied. 

Here, we expect the convergence rate 

t = min{3, — 1 + 2/3 + (1 + a)} = min{3, 2/3 + a}. 

Table [ET4l shows the numerical results for the second order case in which the 
received rates are close to the theoretical expected once. 



a = 2/3 

P 

II 

CO 

a = 1.01 

a = 10/9 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

289 

3.1686e-03 

- 

1.6898e-03 

- 

4.8115e-04 

- 

5.9011e-04 

- 

1.089 

1.2099e-03 

1.39 

6.0844e-04 

1.47 

1.4003e-04 

1.78 

1.5596e-04 

1.92 

4.225 

4.6505e-04 

1.38 

2.1881e-04 

1.48 

3.7277e-05 

1.91 

3.8312e-05 

2.03 

16.641 

1.8057e-04 

1.36 

8.0073e-05 

1.45 

9.9546e-06 

1.90 

9.3965e-06 

2.03 

66.049 

7.0635e-05 

1.35 

2.9545e-05 

1.44 

2.6951e-06 

1.89 

2.3314e-06 

2.01 

263.169 

2.7771e-05 

1.35 

1.0960e-05 

1.43 

7.4481e-07 

1.86 

5.8950e-07 

1.98 

1.050.625 

1.0955e-05 

1.34 

4.0799e-06 

1.43 

2.1075e-07 

1.82 

1.5257e-07 

1.95 



P 

II 

£- 

CO 

a = 3/2 

a = 2.175 

a = 2.275 

DOFs 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

L z -error 

rate 

289 

6.1433e-04 

- 

5.5363e-04 

- 

2.7565e-04 

- 

2.4570e-04 

- 

1.089 

1.5136e-04 

2.02 

1.3540e-04 

2.03 

5.1121e-05 

2.43 

4.1696e-05 

2.56 

4.225 

3.3604e-05 

2.17 

2.8521e-05 

2.25 

7.5320e-06 

2.76 

5.7319e-06 

2.86 

16.641 

7.7002e-06 

2.13 

6.2123e-06 

2.20 

1.1051e-06 

2.77 

7.8407e-07 

2.87 

66.049 

1.8014e-06 

2.10 

1.3642e-06 

2.19 

1.5938e-07 

2.79 

1.0553e-07 

2.89 

263.169 

4.2916e-07 

2.07 

3.0106e-07 

2.18 

2.2723e-08 

2.81 

1.4044e-08 

2.91 

1.050.625 

1.0374e-07 

2.05 

6.6649e-08 

2.18 

3.2138e-09 

2.82 

1.8538e-09 

2.92 



a = 2.375 

DOFs 

L z -error 

rate 

289 

2.2177e-04 

- 

1.089 

3.3912e-05 

2.71 

4.225 

4.3221e-06 

2.97 

16.641 

5.4888e-07 

2.98 

66.049 

6.8762e-08 

3.00 

263.169 

8.5292e-09 

3.01 

1.050.625 

1.0497e-09 

3.02 


Table E.14 L-shaped domain, k = 2: Influence of a for a = 2/3tt and (xo , v/o) = (0,0). 
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E.3 Fichera comer 
E.3.1 Smooth solution 
The geometry is 

n = n F : = (-1, i ) 3 \ [o, i] 3 . 

The discretization is done on lowest order hexahedral elements, regularly refined. 
The exact solution is prescribed to be the smooth solution 

u(x, y, z) = sin((x + y)n) cos(2ttz). 


DOFs 

L2 error 

rate 

316 

0.075444 

— 

3.032 

0.017182 

1.96 

26.416 

0.0039376 

2.04 

220.256 

0.00094597 

2.02 

1.798.336 

0.00023208 

2.01 

14.532.992 

5.7491e-05 

2.00 


Table E.15 Fichera cube. 


E.3.2 Solution of point singularity type 

In the next calculations, the exact solution is given by 

OL 

u = r , 

where r = dist(x,xo) measures the distance from the point xo, which is varied. 
The L 2 -error is computed with a tensor product Gauss rule (5 points in each 
coordinate direction). 



x 0 = (—1, —1, —1), a = 0.55 

x 0 = (0, 0.5,0), a = 0.55 

xo = (0, 0.5,0), a = 2/3 

DOFs 

L2-error 

rate 

L2-error 

rate 

L2-error 

rate 

316 

0.00073994 

— 

0.00069287 

— 

0.00074102 

— 

3.032 

0.00016401 

2.00 

0.00023565 

1.43 

0.00023 

1.55 

26.416 

3.9176e-05 

1.98 

6.8242e-05 

1.72 

6.2591e-05 

1.80 

22.0256 

9.6835e-06 

1.98 

1.9077e-05 

1.80 

1.6589e-05 

1.88 

1.798.336 

2.4305e-06 

1.97 

5.2412e-06 

1.85 

4.3418e-06 

1.92 

14.532.992 

6.1407e-07 

1.98 

1.44225e-06 

1.87 

1.1264e-06 

1.94 


Table E.16 Fichera cube 
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